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Abstract. We describe an elementary algorithm for expressing, as ex- 
plicit formulae in tractor calculus, the conformally invariant GJMS op- 
. . . erators due to C.R. Craham et alia. These differential operators have 

CO I leading part a power of the Laplacian. Conformal tractor calculus is 

I the natural induced bundle calculus associated to the conformal Cartan 

connection. Applications discussed include standard formulae for these 
operators in terms of the Levi-Civita connection and its curvature and 
• a direct definition and formula for T. Branson's so-called Q-curvature 

^ I (which integrates to a global conformal invariant) as well as general- 

Mh , isations of the operators and the Q-curvature. Among examples, the 

operators of order 4, 6 and 8 and the related Q-curvatures are treated 
explicitly. The algorithm exploits the ambient metric construction of 
CO ■ Fefferman and Graham and includes a procedure for converting the am- 

^ ' bient curvature and its covariant derivatives into tractor calculus expres- 

, sions. This is partly based on where the relationship of the normal 

' standard tractor bundle to the ambient construction is described. 

o 
o 
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Conformally invariant differential operators have long been known to play 
an important role in physics and the geometry of many structures related 
I to and including Riemannian and conformal geometries. For example, the 

■ classical field equations describing massless particles, including the Maxwell 

and Dirac (neutrino) equations, depend only on conformal structure |2l [T8]. 
^ I More recently string theory and quantum gravity have motivated several de- 

velopments in mathematics where conformally invariant operators play a key 
role. Many of these could be said to fall under the umbrella of geometric 
' spectral theory where, broadly, one attempts to relate global geometry to 

the spectrum of some natural operators on the manifold. For example, on 
compact manifolds there are programmes to find extremal metrics for func- 
tional determinants of natural operators. Conformally invariant operators 
yield determinants with a workable formula (a so called Polyakov formula) 
for the conformal variation of the determinant thus leading to significant 
progress El El . In another direction there is new progress [35] in relating 
scattering matrices on conformally compact Einstein manifolds with confor- 
mal objects on their boundaries at infinity. This falls within the framework 
of the AdS/CFT correpondence of quantum gravity |iH l IHH j 1^7 1 \'.VA\ . 

In these areas it seems an especially important role is played by natural 
conformally invariant operators with principal part a power of the Laplacian 
A. The earliest known of these is the conformally invariant wave operator 
which was first constructed for the study of massless fields on curved space- 
time. More recently its Riemannian signature variation, usually called the 
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Yamabe operator, has played a large role in the Yamabe problem on compact 
Riemannian manifolds. As an operator on functions it is given by the for- 
mula A — (n — 2)i?/(4(n — 1)), and it governs the transformation of the scalar 
curvature R under conformal rescaling. An operator with principal part 
is due to Paneitz |^ (see also ^BEBl), and then sixth-order analogues were 
constructed in [21 E]- Graham, Jenne, Mason and Sparling (GJMS) solved 
a major existence problem in where they used a formal geometric con- 
struction to show the existence of conformally invariant differential operators 
P2k (to be referred to as the GJMS operators) with principal part A'^. In 
odd dimensions, k is any positive integer, while in dimension n even, k is 
a positive integer no more than n/2. The k = 1 and k = 2 cases recover, 
respectively, the Yamabe and Paneitz operators. 

In dimension 2 the transformation of the scalar curvature can also be de- 
duced from the Yamabe operator by a dimensional continuation argument, 
and the curvature fixing problem corresponding to the Yamabe problem is 
usually known as Gauss curvature prescription. In the late 1980's Branson 
PI irH] observed that the Paneitz operator P is formally self-adjoint and can 
be expressed in the form + ((n — 4)/2)(54, where annihilates constant 
functions and Q4 is a scalar curvature invariant which could play a role par- 
allel to the scalar curvature in higher order analogues of the Gauss curvature 
prescription programme. In dimension 4 the conformal transformation of Q4 
is given by the Paneitz operator, and it follows that the integral of Q4 over 
compact 4-manifolds is a global conformal invariant. On conformally flat 
structures this is a multiple of the Euler characteristic. It has recently been 
established by Graham and Zworski and Fefferman and Graham (H^ ESI [2I| 
that the GJMS operators P2k are formally self-adjoint, and so |3| shows that 
these operators yield an analogous local Riemannian invariant Qn for each 
even-dimensional manifold. There has been considerable recent interest and 
progress in understanding Branson's Q-curvatures, especially in low dimen- 
sions and on conformally flat structures ^| • 

In [22] the GJMS operators are derived from the Laplacian of the ambient 
metric of Fefferman and Graham [23 ESI- This construction is very valuable 
not only in itself but also because of the close links with the Poincare met- 
rics of the conformally compact Einstein theory. On the other hand there is 
another way to generate a conformally invariant operator with principal part 
A'^. The result is usually presented as a simple formulae, first due to M.G. 
Eastwood, as given in (|T5|l . (See [2H| for a derivation and some further re- 
lated developments.) Underlying this formula are two related key tools. The 
first is a geometric construction developed by Eastwood and others [221 
known as the curved translation principle. This construction is a generalised 
and geometric variant of the translation functor due to Zuckerman and oth- 
ers [ISj. The second is a machinery known as tractor calculus [Tl 129 } [TH [T3] . 
This calculus brings the conformally invariant Cartan connection to induced 
bundles and also involves other fundamental conformally invariant operators 
(such as the ones used in this formula). The combination is potent since 
on the one hand it is very easy to expand these tractor formulae in terms 
of the Levi-Civita connection and its curvature (which is useful for the in- 
vestigation of issues such as positivity of the operators), and on the other 
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hand the link with representation theory means one easily obtains rules for 
generalising the operators and how they may be composed with certain other 
conformally invariant operators. See for example (|TH|l . It should be pointed 
out that the tractor formulae are themselves complete and explicit formu- 
lae and can be readily worked with directly without using any knowledge of 
the representation theory aspects. That is essentially the approach below. 
See also [7], for example, where these tractor formulae for conformally in- 
variant powers of the Laplacian are used to construct formally self-adjoint 
conformally invariant boundary problems, higher order conformally invariant 
Dirichlet-to-Neumann operators, and related constructions. 

One problem with the tractor approach up until now has been that, on 
even dimension n manifolds, this had failed to yield the operators of order n 
except for a quotient construction in dimension 4 Here we give a simi- 
lar quotient tractor construction for a sixth-order operator and show that we 
have in fact recovered P4 and Pq. This brings us to one of the main purposes 
of this paper, which is to explicitly relate the tractor calculus approach to 
the GJMS construction. This is achieved in Sectional where an algorithm is 
described for finding a tractor formula for any of the GJMS operators P2fc- 
Remarkably this algorithm does not require solving the Fefferman-Graham 
ambient construction. For low order operators it is essentially trivial and 
quickly recovers the simple tractor formulae for P4 and Pq and yields a cor- 
responding tractor formula for Pg- See Section 14.11 and Proposition In 
Proposition 12.41 we use these formulae to prove directly that these opera- 
tors are formally self-adjoint (verifying directly for these cases the general 
results of (SSIEZI)- Expanding these formulae into formulae in terms of the 
Levi-Civita connection and its curvature simply requires repeated use of the 
Leibniz rule and the definitions of the tractor objects. This is easily auto- 
mated and is done in Section 12.21 The nature of the formulae we use mean 
the calculations have a large number of built-in self-checks which ensure that 
the formulae used are entered and used correctly by the software. Thus over- 
all this demonstrates an effective means to obtain explicit formulae for the 
GJMS operators. It should be pointed out that the formulae in Section [221 
are not in fact the raw output from the expansion of the tractor formulae, 
but rather this output manipulated into the canonical form described in pT] . 
The authors performed these expansions and manipulations mainly by using 
Mathematica and J. Lee's Ricci programme (HHI; this work was performed 
under the assumption of a Riemannian signature metric, but the resulting 
formulae are independent of the signature. 

The most important outcomes of Section 0] are Proposition 14.51 and Theo- 
rem [^J The first of these establishes important features about the form of 
the tractor formulae for the GJMS operators, and the latter exploits this to 
provide some new invariant operators closely linked to the GJMS operators. 
There are several applications of these. One is a direct tractor based con- 
struction of Branson's Q-curvatures. See Proposition 12.71 In fact, this also 
gives a new definition for these invariants. This gives an effective way to cal- 
culate these (Q4 and Qq are treated as examples), and it sheds light on their 
remarkable transformation properties. Another application of Theorem 12.51 
is Gorollarv 12.61 In words this states that except for the k = n/2 case, the 
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theorem yields generalisations of the GJMS operators P2A: that are "strongly 
invariant" in the sense of That is, operators that can be composed with 
tractor bundle valued operators to yield further conformally invariant opera- 
tors. This is one of the key ideas of the curved translation principle. Finally, 
Theorem 12.51 is a crucial ingredient in the general construction in [H] of an 
elliptic conformally invariant operator on 1-forms with close connections to 
the first de Rham cohomology. 

There are other results presented. For example, in Section IT3l we describe 
how to proliferate Riemannian invariants which are not conformally invari- 
ant but have a transformation formula similar to the Branson Q-curvatures. 
These can be viewed as representing terms that could be added to the Q- 
curvature without affecting its key properties and so play a role in generating 
new curvature prescription problems. 

There are also many other potential applications for this work not touched 
upon in this article. For example, the tractor formulae for the GJMS oper- 
ators could immediately be used in a construction parallel to that in [2j to 
produce alternative conformally invariant boundary problems and non-local 
operators based around the GJMS operators. 

It should also be pointed out that the results and ideas in this paper should 
have analogues for CR structures, where one would instead be involved with 
CR-invariant powers of the sub-Laplacian and the ambient construction 
of C. Fefferman [^. The construction presented in this article is in part an 
application of ideas developed in the joint work of one of the authors with A. 
Cap. See ^21 where it is described explicitly how to relate the Cartan/tractor 
approach to the ambient construction of Fefferman and Graham and its 
applications to invariant theory. The relevant aspects of this theory are 
summarised in Section mi There is a corresponding theory for the CR case 

El. 

The authors are indebted to Tom Branson, Audi Cap, Mike Eastwood, and 
Robin Graham for several illuminating conversations. The authors would 
also like to thank the Mathematical Sciences Research Institute and the 
organisers of Spring session in 2001 for helping to make this research possible. 

2. CONFORMAL GEOMETRY AND TRACTOR CALCULUS 

We summarise here an approach to local conformal geometry that is rather 
useful for our applications. This is broadly based on the development pre- 
sented in El) but many of the ideas and tools had their origins in El) H) 
and El- The notation and conventions in general follow the last two sources. 

We shall work on a real conformal n-manifold M, where n > 3. That is, we 
have a pair (M, [g]), where M is a smooth n-manifold and [g\ is a conformal 
equivalence class of metrics of signature (p, g). Two metrics g and g are 
said to be conformally equivalent if 5 is a positive scalar function multiple 
of g. In this case it is convenient to write g = Q^g for some positive smooth 
function Q. Although we assume that the metrics have some fixed signature, 
all considerations below will be signature independent. For a given conformal 
manifold (M, [g]), we shall denote by Q the bundle of metrics. That is, Q 
is a subbundle of S'^T*M with fibre M"*". The points correspond to values of 
metrics in the conformal class. 
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Let 8°" denote the space of smooth sections of the tangent bundle TM, and 
similarly let £a be the smooth sections of the cotangent bundle T*M . In fact, 
we will generally abuse notation and also use these symbols to indicate the 
sheaves of germs of smooth sections and even the bundles themselves. These 
conventions will be carried through to all bundles that we discuss. We write 
£ to denote the trivial bundle over M . Penrose's abstract index notation is 
embraced throughout, so tensor products of these bundles will be indicated 
by adorning the symbol £ with appropriate abstract indices. For example, 
in this notation (S)'^T*M is written £ab- An index which appears twice, once 
raised and once lowered, indicates a contraction. These conventions will be 
extended in an obvious way to the tractor bundles described below. In all 
settings indices may also be "suppressed" (omitted) if superfluous by context. 

The bundle Q is a principal bundle with group M_|_, so there are natural 
line bundles on (M, [g]) induced from the irreducible representations of M+. 
We write £[11:] for the line bundle induced from the representation of weight 
-w/2 on M (that is M+ 9 x x'^'/^ ^ End(M)). Thus a section of £[w] 
corresponds to a real-valued function f on Q with the homogeneity property 
f{x, Cl'^g) = 0"'/(x, g), where is a positive function on M, x E M, and g is 
a metric from the conformal class [g]. We use the notation £a[w] for £a'<^£[w] 
and so on. Note that for consistency with this convention differs in sign 
from the one of [14, Section 4.15]. 

Let £'+[u)] be the fibre subbundle of f corresponding to M+ C M. Choos- 
ing a metric g from the conformal class defines a function / : Q ^ R by 
f{g,x) = where g = 0,'^g, and this clearly defines a smooth section of 
<S-i-[— 2]. Conversely, if / is such a section, then f{g,x)g is constant up the 
fibres of Q and so defines a metric in the conformal class. Thus £"+[—2] is 
canonically isomorphic to Q, and the conformal metric g^^ is the tautologi- 
cal section of £ab\A that represents the map 2] = Q — s- £(^ah)- Prom this 
there is a canonical section g""^ of <S'^''[— 2] such that gab9^'^ = ^a'^ (where 6a'^ is 
the section of £a'^ corresponding to the identity endomorphism of the tangent 
bundle). The conformal metric (and its inverse g"-^) will be used to raise and 
lower indices without further mention. Given a choice of metric g from the 
conformal class, we write Va for the corresponding Levi-Civita connection. 
With these conventions the Laplacian A is given by A = g°'^Vo^b = V^V;, . 
In view of the isomorphism 2] = Q, a choice of metric also trivialises 
the bundles £[w\. In particular we will write for the canonical section of 
£^[1] satisfying g = (C^)~^fir. Conversely a section of £^+[1] clearly determines 
a metric by this relation, so such a is termed a choice of conformal scale. 
This determines a connection on £[w\ via the corresponding trivialisation 
of £[w\ and the exterior derivative on functions. We shall also denote such 
a connection by Vq and refer to it as the Levi-Civita connection. Note in 
particular then that, by definition, Vq^^ = 0, so Va also preserves the con- 
formal metric. The curvature Rab^d of the Levi-Civita connection is known 
as the Riemannian curvature, and is defined by 

This can be decomposed into the totally trace-free Weyl curvature Cabcd and 
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a remaining part described by the symmetric Rho-tensor Pab^ according to 

Rabcd = Cabcd + '^gc[a?b]d + "^Qdlb^a^ci 

where [• • • ] indicates the antisymmetrization over the enclosed indices. The 
Rho-tensor is a trace modification of the Ricci tensor Rob- We write J for 
the trace P"a of P. 

Under a conformal transformation we replace our choice of metric g by 
the metric g = 0,'^g, where is a positive smooth function. The Levi-Civita 
connection then transforms as follows: 

(1) VaUb = VaUb-TaUb-TbUa+gab'^''Uc V aCr = V aCr + wT aCT. 

Here Ub G £b, cr £ £[w], and := ^~^Va^- The Weyl curvature is con- 
formally invariant, that is Cabcd = Cabcd, and the Rho-tensor transforms by 

(2) Pab = Pab - VaTfo + T.Tfe - '^T'^T^Q,,. 

For the density bundle f [1], we have the jet exact sequence at 2-jets, 

where (• • • ) indicates symmetrization over the enclosed indices. Note we have 
a bundle homomorphism £(^ab) [1] ~^ given by complete contraction with 

gab 'pj^jg jg gpin via p I— > ^pQab ^nd so the conformal structure decomposes 
<5(Qfe)[l] into the direct sum <S((jfe)p [!]©£'[— 1]. Clearly then iS(afe)o [1] is a smooth 
subbundle of J^(f [1]), and we define £^ to be the quotient bundle. That is, 
the standard tractor bundle £^ is defined by the exact sequence 

(3) 0^f(,,)Jl]^ J2(f[l])^^^^0. 

The jet exact sequence at 2-jets and the corresponding sequence at 1-jets, 
viz — > £^a[l] — > •^^(^[1]) '^'[1] ~^ 0, determine a composition series for 
£'^ which we can summarise via a self-explanatory semi-direct sum nota- 
tion £'^ = f [1] 6-fa[l] 1]. We denote by the canonical section of 
£^\V] := £^ ® £[\] corresponding to the mapping <?[—!] — > f 

Composing the canonical projection J^(f [1]) £^ with the 2-jet operator 
p yields an invariant differential operator : £[1] ^ £^ . On the other 

hand, if we choose a metric g from the conformal class, then the map 

jla ^ [^^DMx)]a ■■= (^(x), Vaa{x), -i(A + J)a(x)), 

induces an isomorphism £^ — > <S[l]©<Sa[l]®<?[— 1] =: [£^]g of vector bundles. 
Tautologically the displayed formula for ^[D^a{x)]g gives the operator 
in terms of this decomposition. If the image of G £^ is [V^]g = {a, lJ,a,T), 
then from the change in the Levi-Civita connection Q we get 

[V\ = (a^r) = {a, fia + crTa, r - Tbfi' - iaT^T^). 

This transformation formula characterises sections of £^ in terms of triples 
in £[1] £a[l] (B£[-1]. With a fixed rescaling of the map £[-1] ^ £^, we 
have [^"^Ig = (0,0, 1). It is convenient to introduce scale-dependent sections 
^Ab g £^^''[—1] and G £""^[—1] mapping into the other slots of these 
triples so that \y^]g = {cr,fJ,a,'T) is equivalent to 
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Figure 1. Tractor inner product 

If and are the corresponding quantities in terms of the metric g = 
0^(7 then we have 

(4) Z^^ = Z^^ + T*X^, Y^ = Y^- TbZ^^ - ^TbT^X^. 

The standard tractor bundle has an invariant metric Hab of signature 
{p + l,q + 1) and an invariant connection, which we shall also denote by 
Va, preserving Hab ■ If is as above and S £^ is given by [y_^]g = 
(£,/fb,z), then 

hABV^Vf = fi^fi^ + err + Ta. 
Using hAB and its inverse to raise and lower indices, we immediately see that 

YaX^ = 1, ZAbZ'^c = Qbc 

and that all other quadratic combinations that contract the tractor index 
vanish. This is summarised in Figure ^ Thus we also have YaV"^ = 
T, XaV^ = a, ZAbV^ = and the metric may be decomposed into 
a sum of projections, Hab = Za'^Zbc + XaYb + YaXb ■ 

If for a metric g from the conformal class V"^ € 6"^ is given by [V^]g = 
{a,fia,T), then the invariant connection is given by 

/ VaCT - /ia \ 

(5) [VaF^lg = ^af^b + QabT + PabC^ . 

The tractor metric will be used to raise and lower indices without further 
comment. We shall use either "horizontal" (as in [V^]g = {a,fib,T)) or 
"vertical" (as in notation, depending on which is clearer in each given 
situation. 

Tensor products of the standard tractor bundle, skew or symmetric parts 
of these and so forth are all termed tractor bundles. The bundle tensor 
product of such a bundle with £[w], for some real number weight w, is termed 
a weighted tractor bundle. For example £aiA2-Ai [w] = £ai 'Si • • ■'Si£ai'Si£[w] 
is a weighted tractor bundle. Given a choice of conformal scale we have the 
corresponding Levi-Civita connection on tensor and density bundles. In 
this setting we can use the coupled Levi-Civita tractor connection to act on 
sections of the tensor product of a tensor bundle with a tractor bundle. This 
is defined by the Leibniz rule in the usual way. For example if u^'V'^'a G 
£''^£^ S £[w] =: £''^[w], then VaU^V^a = {VaU^)V^(J + u\VaV^)(J + 
u'^V^'VaO'- Here V means the Levi-Civita connection on G £^ and a G 
£[11}], while it denotes the tractor connection on V*-^ G £'-^ . In particular with 
this convention we have 



(6) VaXA = ZAa, V aZAb = -PabXA - YAOab , ^ aYA = abZ , 
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which for the purposes of automating calculations is a very useful description 
of the tractor connection. 

Note that if F is a section of <f*[u;], which means simply some tractor 
bundle of weight w, then the coupled Levi-Civita tractor connection is not 
confomally invariant but transforms just as the Levi-Civita connection trans- 
forms on densities of the same weight. That is 

under the conformal rescaling 

g ^ g = (cf. ([TJl ) . It is an elementary 

exercise using the last transformation formulae and Q to show that, for 
F G <S* [w] , the formula 

(7) D^^V := 2wX^Py^W + IX^^Z^^^V^V 

determines an invariant operator D^^ : £^[w\ —>■ £^[^^1 <f^*[iL']. (This 
was first developed in early versions of [22] and is closely related to the 
"fundamental D" operator developed in [l^-) Since we can vary the weight 
and the tractor bundle £^ , D^^ is really an entire family of operators. The 
point is that with the way we have defined V, the same formula works for 
the entire family, and so it is reasonable to let the single symbol D^^ denote 
all of these operators. We abuse terminology and describe it as an operator. 
(The Levi-Civita connection is usually used this way.) If we have a single 
formula Op that gives a family of conformally invariant operators 

Op : ^* (g) V ^ £:* ^ w 

as we range over all tractor bundles £^ then, following P^], we describe Op 
as a strongly invariant operator. For example D"^^ is strongly invariant. 
As already pointed out D^^ is rather more universal since we can vary the 
weight w as well. Thus we can form compositions of this operator with 
itself, and in particular we consider h^^ Dj^(^QD\B\py^V for V G some 
weighted tractor. Expanding this out using 0; and the Leibniz rule for 
V, it is easily verified that it may be re-expressed in the form 

h^''Dj,^QD\s\P),V = -X^qDp^^V, 

where D is some operator determined explicitly in the calculation. Since the 
map £p[w — 1] — > £i^pQ)g[w] given by Sp i-^ Xi^qSp^^ is injective, this estab- 
hshes Da : £'^[w] ^ £a®£'^[w- 1] conformally invariant differential 
operator on weighted tractor bundles. For V G ^"^[tt;], this is given by 

(8) D^V := {n + 2w- 2)wY^V + {n + 2w- 2)Z'^^VaV - X^DV, 
where 

(9) nV := VpVPy + wW. 

So Da is in fact precisely the tractor D-operator in Note the identity 

(10) DaX'^V = in + 2w + 2){n + w)V, 

which we will use later. 

The curvature O of the tractor connection is defined by 
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for F*^ G £'~^ and is precisely the local obstruction to conformal flatness. 
(That is locally there is a flat metric in the conformal class if and only if this 
curvature vanishes.) Using (EJ and the usual formulae for the curvature of 
the Levi-Civita connection we calculate (cf. ^) 

(11) ^abCE = Zc'^ZE'^Cabce — ^^^IC^ E]'^ [aPb]e- 

It is straightforward to use this and © to show that if F € £ce---f[w\., then 

[Da,Db]Vce-f = 
^^2^ {n + 2w- 2) [Wabc^Vqe...f + 2wnABC^VQE-F 

+AX[A^Bfc'^'^sVQE-F + ■■■ + Wabf'^VcE-Q 
+2wnABF^VcE-Q + 4X[A^BfF'^^sVcE-Q]. 

Here 0.abce = Za"" ZB^^ahCE, ^BsCE = ZB^^bsCE, and 

(13) WabCE = (n - 4:)nABCE - 2X^AZB\^''^pbCE- 

It follows that on conformally flat structures \D a^Gb\^ce---f = 0. Similarly 
it is easily verified that \Db^Gc\ annihilates densities. 

We should point out some features of Wabce- Firstly, it is conformally 
invariant. One can already see this from ifT^ by setting w = and then con- 
sidering sections Vc of £c such that VaVc vanishes at a given point. This 
is also immediately clear from the formula Wab^ L '■= X^pQab]^ L 

(see [Sni), which is readily verified. From this several things are immedi- 
ately clear. Firstly, Wabce vanishes on conformally fiat structures. Next, 
we have that Wabce = W^[yi_B][C£;] and that it is trace-free ( since ^abce is 
annihilated by contraction with X^ on any index). Furthermore expanding 
lfT3)l reveals that W[abc]e = 0- Thus Wabce has "Weyl tensor symme- 
tries". Whence it is immediately clear that Wabce is also annihilated upon 
contraction with X^ . 

Finally we should comment on the uniqueness of this tractor calculus. In 
sections 2.6 and 2.7 of ^3] it is shown that the transformation properties (@J 
and the form of the connection ^ identify £^ and its tractor connection Va 
as above as a normal tractor bundle and connection corresponding to the 
defining representation of SO(p-|- l,q + 1). Let V be IR"'^^ as the representa- 
tion space for the standard (or defining) representation of SO(p-|- l,q + 1). 
We can construct ^1] from the pair (<S^, Vq) a principal bundle Q which is 
the frame bundle for f corresponding to the metric and filtration. This 
has fibre P, a certain parabolic subgroup of SO{p -|- 1, g -|- 1). A Cartan con- 
nection u; on ^ is determined by V. This is the normal Cartan connection 
on Q such that is the vector bundle connection induced from uj. That is 
the normality condition on the pair {£^,'S/a) is equivalent to the pair {G,uj) 
being a normal Cartan bundle and connection in the sense of [15]. 

2.1. Conformally invariant powers of the Laplacian. Since the tractor- 
D operator constructed above is well-defined on any weighted tractor bundle, 
we can compose the tractor-D operators. It is clear from the formula for the 
tractor-D operator that any such composition will yield a natural operator. 
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that is an operator which can be written as a polynomial formula in terms of 
a representative metric, its inverse, the metric connection and its curvature. 
On densities of the appropriate weight and with some minor adjustment a 
composition of this form will lead to conformally invariant operators with 
principal part a power of the Laplacian. 

First let us observe how the conformal Laplacian arises from the tractor 
machinery. Let / G £[1 — n/2]. Then observe that immediately from (jH)) we 
have DaI = —XaOJ- Since Da is conformally invariant we have immedi- 
ately that for / G f [1 — n/2], □/ is conformally invariant. From (jHl this is 
(V^Va + ^^J)/ - the usual conformal Laplacian. 

Now suppose that instead we have / G f *[1 — n/2], where here and below 
£^[w] will be used to indicate any tractor bundle of weight w. We still have 

(14) DaI = -XaOJ, 

but now in □/ = (V^Vq + J)/, V means the Levi-Civita tractor coupled 
connection. In particular this establishes a strongly invariant generalisation 
of the Yamabe operator on tractor sections of the said weight. 

It is clear from our observations that that there is a conformally invariant 
operator 

ODa,Da,---Da,_, : £[k - n/2] ^ Sa,A,-A,_A-1 - n/2]. 
In the conformally flat case this already yields an operator between densities 

(cf. EHl). 

Proposition 2.1. On conformally flat structures, if f & £[k — n/2], then 

nDA,_, ■ ■ ■ DaJ = {-i)''''Xa, ■ ■ ■ XA,_,a2kf, 

where \I\2k '■ £[k — n/2] £[—k — n/2] is a conformally invariant operator. 
Locally we can choose a flat metric from the conformal class. This determines 
a connection in terms of which we have 02kf = A^/. 

Proof. In any choice of conformal scale, expand out \I\Da^_-^ ■■■DAif via 
the formula ^ and move the X,Y,Z^s to the left of all V's via the iden- 
tities (inj. It is immediately clear that the highest order term is precisely 
{—1)''~^Xai ■ ■ ■ XAf,_i^'^f and that any other coefficient of Xa^ ■ ■ ■ -^Afe_i 
involves the curvature Pab or its trace. 

On the other hand, on conformally flat structures, [Da, Db]V = for 
V any weighted tractor field. Thus Da^ • ■ ■ DA^f is completely symmetric 
for any I G In particular ^Da^_^ - --DAif G £(Ai---Ak-i)[~'^ ~ n/2] 

and Da^Da^_^ ■■■DAif G £(^Ai---Ai.)[~n/2]. Consequently it must be that 
= D[A^:DA^:_-^^ ' ' ' DaJ- But • • • DaJ' has weight 1 - n/2, so from 

lfn)l this implies Xj^^DDy^^^jD^j. ^ •••Dai/ = 0. It follows immediately 

that nDAk_^DAk_2 ■ ■ ■ DaJ = i-l)''~^XAi ■ ■ ■ XA,^_-,02kf for some opera- 
tor □2A:- With the above we are done. □ 

If we are happy to work in the scale of a flat metric then there is an even 
simpler proof along the lines of the proof of Proposition 14., SI We leave this 
for the reader. 
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By the same ideas as in the proof above, it is easy to use ifT^ and ^ to 
show that, if A; > 3, then X[^^nD^^_^]Z?Afe_2 ■ • • -Dai / / for / G £[k-n/2] 
on a general conformally curved manifold. Thus the proposition fails if we 
remove the requirement of conformal flatness. One way to generalise the □2^ 
is as follows. 

Consider 

for / € f *[A; — n/2]. This is manifestly strongly conformally invariant in all 
dimensions and for all positive integers k. Furthermore by the identity IjlOjl 
we have that, on conformally flat structures, 

(15) • • • D^'^-^DDa,., ■■■DaJ= (j[{n - 2i)ii - 1) j D^kf. 

On the other hand for general conformally curved structures, suppose that 
the dimension n is odd or satisfies 2k < n. Then we can define ^2kf by (|TK|l . 
and this gives a conformally invariant operator 

(16) Dafc :^*[A;-n/2] ^^*[-A:-n/2] 

with principal part A'^. Here, as usual, £'^[k — n/2] indicates any tractor 
bundle of weight k — n/2. In these dimensions this generalises the operator 
□2fc of the proposition. Although we do not wish to describe the curved 
translation principle (221 El, it is worth pointing out that it is partly illus- 
trated here. The tractor formula lfTK]l manifests as a "translate" of the 
Yamabe operator □. In fact proceeding in smaller steps it demonstrates ^2^ 
as a translate of 02^-2 • 

Before we move on, let us demonstrate that the operators 02^ are formally 
self-adjoint. We summarise some results we need from Section 7 of [2j in 
the following proposition. These results can be verified easily using the 
definitions above, and it is important for our needs to note that this works 
in a rather formal manner. That is, we can leave the dimension and weight 
as unknown in the calculations. 

Proposition 2.2. On a conformal manifold M we have 

(i) If ijj^ € r£''^[?i;] and ip € r£^[l — n — w\ is compactly supported on M , 
then 

[ ipDAiJ^= [ {DA^)i''^. 
Jm Jm 

(ii) If £^ is any tractor bundle, then f* is canonically isomorphic to its 
dual £^ via the tractor metric and for any pair ip^,ip^ G r£'*[l — n/2], 
(£^[1 - n/2] := £^ £[1 - n/2]) we have 

Jm J 

Since n] is naturally identified with the space of volume densities the 
integrals are well-defined. Now part (ii) of the proposition asserts ^2 = □ is 
formally self-adjoint, while the same result for D^^ ■ ■ ■ D'^''-^UDa^,_^ ■ ■ ■ Da, 
follows immediately from this and repeated use of (i). So the formal self- 
adjoint property of □2^ is proved. It should be pointed out that as well 
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as observing that D^^ ■ ■ ■ D'^''-'^\I\Da^._i ■ ■ ■ Dai/ recovers a conformally in- 
variant power of the Laplacian, M.G. Eastwood also observed the formal 
self-adjoint property. It is clear from ifT^ that, unfortunately, this formula 
does not yield a conformally invariant operator of order n on even dimen- 
sional structures, yet the existence of such an operator is guaranteed by the 
construction of |32| . 

Recall from Section |21 that if / G f'Mi then [DA,DB]f = 0, and so the 
k = 2 case of the proposition does hold on general conformal structures. In 
particular, as observed in |2Bj, for / G £"[2 — n/2] we can define P4/ by the 
quotient formula 

BDAf = -XaP'J- 

Then P'^ has principle part A^. This construction works even when n = 4, 
and in other dimensions P4 = ^4 as defined above. It is not hard to do the 
next even order in a similar way. If now / S <5[3 — n/2], then [Db, Dc]f = 
and hence DsDcf = Df^B^c)! ■ Now it is a short exercise, using lfT2|l and 
the definition of Wabcd once more, to show that 

(n - A)X[AnDB]Dcf = -2X[AWB]^jDsDTf. 

Now from the fact that DsDxf is symmetric and that Wbsct has Weyl 
tensor type symmetries, we can deduce that in dimensions n 7^ 4, 

(17) PBcf ■■= DDBDcf + -^WB^c'^DsDrf 

n — 4 

is symmetric (i.e. PbcI € £(^bc)[~^ ~ ^/'^])- On the other hand, from the 
previous display -Y[a-Pb]c/ = 0. Thus, for n / 4, 

PBcf = XbXcPU, 

where Pg is a conformally invariant operator f [3 — n/2] £[—3 — n/2] 
generalising (for the allowed dimensions) the sixth-order operator of Propo- 
sition In particular this works in dimension 6. 

We should point out that although \I\2k as defined by (|T^ is manifestly 
strongly invariant, we cannot conclude this for Pg as defined here. The 
operator Pbc defined above is clearly invariant when acting on weighted 
tractors, but the argument here to deduce that PbcI has the form XbXcPqJ 
relies on the vanishing of [Da^Db]! ■ 

We will establish in Section (see in particular Subsection I4.1|l the fol- 
lowing result. 

Proposition 2.3. The operators P^ and Pq defined by the tractor expressions 
above are precisely the fourth-order and sixth-order GJMS operators. That is 
P4 = P4, Pq = Pq. a tractor expression for the eighth-order GJMS operator 
Ps is as follows: 

XAXBXcPgf = 

-UDADBDcf - ^^WA^'B^DpDQDcf - ^^Wa"" DpDsDQf 
-^^XAUB^'c'^DpDQf + ^-^^XAD^WB^'c'^DEDpDQf 

+4 („_4V„-g) XAWB''c'^Wp'Q''DsDTf, 
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where all operators act on everything to their right, in a given term, and 
Ub^C^ is the tractor field 

^ {W^Pb^Wfac^^ + W^Pc^Wbaf"^ + W^pQ^Wbacf) ■ 

Here D^Wb^ c'^ DeDpDqJ is taken to mean {Wb^ DeDpDqJ) (and 
not {D^Wb^ c^)DEDpDQf). This and similar conventions for other oper- 
ators and situations will apply throughout the paper. 

We should emphasise at this point that the tractor formulae for P4, Pg 
and Pg above, and similar ones for the higher order P2k that we could easily 
construct via the algorithm of Sectional are genuine formulae for the GJMS 
operators. No further algorithm is required. They are valid on any confor- 
mal manifold where the given GJMS operators exist. In this tractor form 
they are already suitable for many applications, such as establishing strong 
invariance or constructing related operators. The remainder of the section 
will demonstrate this. 

We begin by using the tractor formulae directly to show that the operators 
Pi-,P& and Ps are formally self-adjoint (FSA). We treat these in order. For 
f e£[2- n/2], we have 

(18) DDAf = -XaPaI which implies D^DDaI = (n - 4)Pif. 

We have already observed that D^^ ■ ■ ■ D'^''-'^ODA^^_^ ■ ■ ■ Da^ is FSA on 
£[k — n/2]. So from the second of these it is clear that P4 is FSA in dimen- 
sions other than 4. From the expressions IHI) and © it follows that D^DDAf 
and nDAf, as expressions in terms of Levi-Civita covariant derivatives of /, 
Pab and J, are polynomial in n. So from (|TH|l it is clear that (4 — n) divides 
this expression for D'^ODAf and so P4/ is also given as a formula polyno- 
mial in n and the Levi-Civita covariant derivatives of /, Pab and J. Working 
among tensors of this form a calculation to verify the FSA property of P4 
(in dimensions greater than 4) can be carried out formally in dimension n, 
since Proposition 12.21 is established that way. It follows immediately that 
the same calculation must work when we set n = 4. Thus P4 is also FSA 
in dimension 4. Now for Pg, let / G £^[3 — n/2] and note that □L's-Dc/ 
and Wb^ DsDxf are polynomial in n. Thus PbcI is rational in n with 
a singularity only at n = 4. From PbcI = ^sXcPef we have 

(n - A)D^D''PBcf 

= {n- A)D^D^DDBDcf + 2D^ D^Wb^ DsDrf 

= 2(n-4)2(n-6)P6/. 

Now since Wbsct has the Weyl tensor symmetries (in fact here we just 
need Wbsct = Wctbs = Wtcsb), it follows from Proposition 12.21 that 
D^D^Wb^c^DsDt is FSA on £[3 - n/2]. We know D^D^UDbDcI is 
also FSA and as expressions in terms of Levi-Civita covariant derivatives 
of /, Cabcd, Pab, and J, both of these and {n — 4)Pbc are polynomial in 
n. Thus the expression like this for Pbc f is divisible by (n — 6), so 

reasoning as for P4, we quickly conclude that Pg is FSA in all dimensions for 
which it is defined. Finally, since Uabcd also has Weyl tensor symmetry (as 
readily verified directly or since it corresponds to A.Rabcd as in Sectional), 
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it follows that D^D^Ub^c^DpDqJ is FSA for f e £[4: - n/2]. A similar 
comment applies to the other terms on the right-hand side of the formula, 

6(n - 4)2(n - 6)(n - 8)^3/ = 

(n - A)D(^D^D^nDADBDcf + 2D^ D^Wa^ DpUgDcf 

+2D^D^D^WA''c'^DpDBDQf - 4^J^D^ D^Ub"" c"^ DpDq/ 

+2^D^D^D^WB''c^DEDpDQf 

^A^D^D^WB^'c'^Wp^Q^DsDTf, 

which follows from the earlier display for P%. This shows immediately that 
P%f is FSA in dimensions other than 8, and then, arguing as in the previous 
cases, we can deduce that it is also FSA in dimension 8. We have directly 
proved the following. 

Proposition 2.4. The GJMS operators -Rjj-Pe o,nd Pg are formally self- 
adjoint. 

In fact the result is also immediate from the formulae for these operators 
in Section 12.21 The formulae there are given in terms of the Levi-Civita 
connection and its curvature and are in a canonical form that manifests the 
formal self-adjoint symmetry. (It should be pointed out that in deriving 
those formulae formal self-adjointness was not assumed.) 

Recently the entire family of operators P2k have been shown to be formally 
self-adjoint by other means p5l i27i | . 

In Section we will show that there are similar tractor formulae for all 
of the GJMS operators and that these tractor formulae share some of the 
qualitative features of the examples above. In particular we will prove the 
following theorem. 

Theorem 2.5. (i) In each dimension n and for each integer 2 < k, with 
k < n/2 if n is even, there is a conformally invariant differential operator 

EcD^^ ■.£AB[k-2- n/2] ^£cD[2-k- n/2] 

such that 

Ecn'^^'DADBf = XcXDP2kf 
for f G £[k — n/2], where P2k '■ £[k — n/2] £[—k — n/2] is the order 2k 
GJMS operator. The operator is given by a formula which is a partial con- 
traction polynomial in Da, Wabcd, Xa, hAB and its inverse 
(a) With n and k as in (i) and for f* any tractor bundle, there is a confor- 
mally invariant differential operator 

EcD^^ : £'^ (S) £AB[k - 2 - n/2] ^ £'^ £cd[2 - k - n/2] 

which generalises the operator of part (i). 

(Hi) On conformally flat structures the operator E is, up to a non-zero scale, 
□2fc-4- In this case, given a choice of fl,at metric from the conformal class, 
E is, up to a non-zero scale, A*^"^. 

Proof. We have already observed that for / E £^[2 — n/2], we have (see (CHI)) 
UDAf = -XAPif. Since then G £a['\- - n/2] we have XBUDAf = 

—DsDAf, and so we have (i) for k = 2. Otherwise establishing part (i) is 
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the primary purpose of Section ^ More precisely, from the discussion there 
we obtain Proposition I4.5( which asserts that 

Xa,_, ■ ■ ■ XA,P2kf = {-lf-^UDA,_, ■■■DaJ + ^A,.r-A/'^DpDQf. 

Applying D^^ ■ ■ ■ D^fe-i to both sides of this and using (fT?1|l . we obtain 

Xa.XaAI = EA,A/^DpDQf, 

where 
(19) 

Ea.a/"^ = U^(k-2)^Af6A,^ + 

As explained in Proposition 14.51 ^ is given explicitly by a sum of terms 
each of which is a monomial in Da, Wabcd, Xa, hAB, and its inverse h^^ . 
Each such monomial is thus a composition of strongly conformally invariant 
operators. So Ea^A2^'^ is a sum of compositions of strongly conformally 
invariant operators. Just knowing that Ea^A2^'^ is a sum of compositions 
of conformally invariant operators of this form gives part (i). Then part (ii) 
is immediate from the fact that these are strongly invariant operators. 

Next we show part (iii). Prom Proposition 14. 51 each term in the expression 
for "if is of degree at least 1 in Wabcd- The latter vanishes on conformally 
flat structures. Thus, from (|T9|l . on such structures E is just □2^-4. From 
Proposition I2.H given a choice of flat metric from the conformal class, we 
have □2fc-4 = A''-^. □ 

Remarks: In regard to part (i) of the theorem we should point out that 
Section 0] not only establishes the existence of a formula for E which is a 
partial contraction polynomial in Da, Wabcd, Xa, hAB, and its inverse 
/i^^, but describes an algorithm for finding such a formula. 

It seems likely that the operator E in the theorem is formally self-adjoint. 
Note for example, if we write E* for the formal adjoint of E, then from 
Proposition 12.21 the identity lfT?Hl . and the result that P2A: is formally self- 
adjoint, we have iUi=k-iin-2i)ii-l))P2k = D^D^ E^^^Q DpDg on £[k- 
n/2]. 

Finally we should add that the terms which distinguish the P2A: from the 
□2fc do not vanish in general. At least we have verified by direct calculation 
that for / G £[3-n/2] (and n 7^ 4) the leading term oiD^ D^Wb^ c'^ DsDrf 
is a non-zero scalar multiple of C"" cdeC^^'^'^'^ bf ■ Thus Pq, is not simply a 
scalar multiple of Dg. 

It is a non-trivial matter to know when conformally invariant operators 
have strongly invariant generalisations. Some do not. For example in dimen- 
sion 4 we know there is a conformally invariant operator P4 : £" — £"[—4] with 
principal part A^. Suppose there were a strongly invariant generalisation of 
this. Then, in particular, it would give a conformally invariant operator 
Ha^ : £b — 4] with principal part A^. (Here we mean the principal 
part as an operator between the reducible bundles indicated.) Then, in the 
case of Riemannian signature conformal 4-manifolds, using the ellipticity of 
this, (fTT)]) . Proposition 12.21 and the differential operator existence results in 
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the conformally flat setting (cf. |^) we can conclude that we would have 
a conformally invariant operator D^Ha^Db ■ £[1] ^[—5] with principal 
part (on arbitrary conformal 4-manifolds). This contradicts C.R. Gra- 
ham's non-existence result and so we can conclude the operator Ha^ 
does not exist. (See also j21]-) However P4 does have a strongly invariant 
generalisation in all other dimensions. This is just ^4 as a special case of 
lfT?)]l . More generally, a consequence of part (ii) of the theorem is that the 
GJMS operators admit strongly invariant generalisations except in the 
critical dimension n = 2k. That is, we have the following proposition on 
n-dimensional conformal manifolds: 

Corollary 2.6. For each integer k > 1, with 2k < n if n is even, there is 
a (tractor) formula that gives, for each tractor bundle 8^ , a formally self- 
adjoint differential operator 

Pl:£%-n/2]^8^[-k-n/2l 

where £'^[w\ := ® £[w]. The operator has principal part A'^ and can he 
expressed as a sum of '0\2k and a contraction polynomial in Da, Wabcd, 
Xa, hAB, (ind its inverse h^^ . In the conformally flat case the operator is 
□2fc • In the case that then = P2A; • 

Proof. Since Da is strongly invariant and since also, from (ii) of the theorem, 
Eab^^ is strongly invariant, it follows that there is a conformally invariant 
operator (F := (ntfc-i(»^ " 2i)(i - 1))-^D^D^Eab'''^DpDq) : - 
n/2] — > £■*[— A: — n/2] for any tractor bundle. By ifTTHl this precisely recovers 
the GJMS operator P2k if £'^[k — n/2] is simply the density bundle £[k — 
n/2]. Now consider the formal adjoint F* of F. This is another conformally 
invariant operator £^[k — n/2] £^[—k — n/2] (where as usual we identify 
with its dual via the tractor metric). Since P2k is formally self-adjoint, 
it is clear that, when applied to f [/c — n/2], F* also recovers the GJMS 
operator. Thus {F + F*)/2 is the required formally self-adjoint operator. 

It is clear from Proposition 14.51 that we can express F by a formula which 
is a sum of □2fc and a contraction polynomial in Da, Wabcd, ^a, hAB, 
and its inverse h^^ . From that proposition we also have that each term in 
the latter polynomial expression is of degree at least 1 in Wabcd- Using 
Proposition 12.21 and the formal self-adjoint property of \^2k, we see that 
there is an expression for F* as a sum of and a contraction polynomial 
in Da, Wabcd, Xa, hAB, and its inverse h^^ . Again each term in the 
latter polynomial is of degree at least 1 in Wabcd- So the final part of the 
corollary follows from these observations and Proposition 12.11 □ 

2.2. Conventional formulae. There are circumstances where it is useful 
to have explicit formulae for the GJMS-related operators and invariants in 
terms of the Levi-Civita connection and its curvature. These formulae are 
generally cumbersome. But the various curvature terms are closely related 
to the spectrum of the operator, so it is important to be able to extract these 
explicitly. In particular, for example, issues of positivity can be investigated 
directly in this setting. Moreover such formulae are ready to be mechanically 
rewritten in local coordinates should this be required. 
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Here we will describe how to re-express tractor formulae for P2kf into 
formulae which are polynomial in g, its inverse, V (meaning the Levi-Civita 
connection), C, P and J, and of course linear in /. 

For the most part, the process is simply an expansion of the tractor for- 
mulae using the definitions above. Consider the Paneitz operator P4 first. 
We observed in Proposition O that for f € £[2 - n/2], -XaP^ = UDAf. 
So P4/ = —V^UDAf, and we could simply calculate this scalar quantity 
Y'^\JiDAf ■ In fact we prefer to expand the entire tractor valued quantity 
□D^/ using ((HI) and ©. According to its definition, Da lowers weight by 
1. Thus D-Da/ is given by 

(VftV" + (1 - ^)J)((4 - n)YAf + 2Za"V„/ - Xa{V,V' + (2 - ^)J)/). 

Now we simply move the Xa-,Ya and Za"' to the left of all other operators 
by repeated use of This is easily done by hand and simplified via the 
Bianchi identity to yield 

-UDaI = XA(A2/-(n-2)JA/ + 4P^^V,V,/-(n-6)(VM)VJ 

- ^(AJ)/ + M!^j2/ _ _ 4)P,,P^V). 

The coefficient of Xa on the right-hand side is a formula for the Paneitz 
operator. Note that the coefficient of Ya and the coefficient of Z^" both 
turned out to be zero. Of course this is exactly as predicted by our formula 
—XAPif = D-Da/, but it provides a very useful check of the formulae to 
verify this. So this is all there is to producing the required formula for i-*4 
from the tractor formula. Before we continue with the general case let us 
just reorganise the result. 

For any linear differential operator on densities of the appropriate weight 
there is a canonical form for the formula which, among other features, mani- 
fests the symmetry in the formally self-adjoint and formally anti-self-adjoint 
parts |l2T]. As already observed, the Paneitz operator is formally self-adjoint. 
Applying this idea to the formula above yields 

(20) P4/ = V,V,Sf' V,Vz/ + vSi^.f + ^QiJ. 

Here S^'"' and S^^' are the tensors {l/?,){g^^g^^ + g^'^g^^ + g'^g''^) and 




respectively, and Q\ „ denotes the scalar 

(21) §j2-2PijP^^-AJ. 

In (pn|l . the V's act on all tensors to their right within the given term. 

We now discuss the general case. Explicit tractor formulae are readily 
produced by the algorithm described in Section^ and so we shall suppose 
we are beginning with a formula for P2A: as described in Proposition EHl The 
formulae for and above (see Proposition I2.,S|1 give explicit examples 
that can be kept in mind. These formulae are polynomial in □, Da, Wabcd, 
Xa, hAB, and its inverse h^^ . We replace each of these with its formula 
in terms of the coupled tractor-Levi-Civita connection V, Xa and so forth 
according to the formulae (jHl, lfTI|l - and l[T^ . In doing so we note that 
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W has weight —2 and that D lowers the weight of a tractor by 1. Next we 
move all occurrences of and Ya to the left of the V via repeated 

use of ((HI). At the end of this process all tractor valued objects are to the 
left of the remaining V's, and so at this point these V's are simply Levi- 
Civita covariant derivative operators. Next we use the inner product rules of 
Figure n to simplify the resulting expression. The formula for P2kf is then 
simply the overall coefficient of XaiXa2 ■ • Prom Proposition 14.51 all 

other slots of the tractor expression vanish. That is, the sum of the terms 
that do not contain XaiXa2 • • • XAk_i is zero. Verifying this or even partly 
verifying this provides a very serious check of all formulae and any software 
that are used in the calculation. For example, one can verify that the sum 
of the terms containing Zai'^Xa2 ■ ■ ■ ^Ak_i vanishes. 

This procedure is very simple. But there are many terms involved, as the 
next examples will illustrate. Thus it becomes very useful to be able calculate 
via a suitable computer algebra system. For the examples below the authors 
used Mathematica and J. Lee's Ricci program which proved to be very 
effective. Certainly the Pg case is beyond a reasonable hand calculation. 
The use of software and the self-checking nature of the formulae as discussed 
above mean that one can be confident of the final result. 

As a technical point for these calculations, we describe a simple technique 
which can considerably reduce the computing time they require. One can 
implement this technique by developing a short computer programme. We 
begin by noting that certain steps in the computation may produce tractor 
inner products of the form BiB2...Bi^^^^^"'^^ , where the indices Bi, B2,. . . , 
and i?£ appear as subscripts or superscripts attached to the tractors y, Z, 
and X. Suppose that £ is large and that the tractors ^BiB2...Bi and $^1^2. ..-B^ 
are the sums of many terms. Suppose also that no derivatives of Z, or 
X occur. The tractor "^BiB2...Bi is a linear combination of the following 3^ 
terms: 



yB,YB2- 


■Yb,_ 


.Yb, 


Yb,Yb2 ■ 


■Yb,_ 




Yb,Yb2- 


■Yb,_ 


i^Bt 


Yb,Yb2- 


■ Zbi,_ 


-1 ^ Bi 


Yb,Yb2- 


■ ^Be. 




Yb,Yb2- 


■ Zbi,_ 


^'^-^Xn 
-1 ^^Bi 



XbiXb2 ■ ■ ■ Xb^_^Xb^ 

The coefficients of these terms may, of course, be very complicated. By 
raising indices we may write ^^i^z.-Bf similar linear combination. Each 
term of each linear combination may be paired off with at most one term 
in the other linear combination so as to give a nonzero inner product. We 
compute the 3^ possible inner products and add the results. 

We conclude this section with the calculation of P% and Pg via these 
methods, beginning with the tractor formulae indicated in Proposition 12.31 
As a check, the authors verified the vanishing of the overall coefficient of the 
Zb^Xc term in the expansion of lfT7)l . In a similar fashion, they also verified 
the vanishing of the overall coefficients of the XbZc^ and Za^XbXq terms 
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in the expansions for Pq and Ps, respectively. This involved the use of the 
Bianchi identities, tensor symmetries, and changes in the order in which 
covariant derivatives are taken. The authors also manipulated the resulting 
formulae for the GJMS operators into the canonical form suggested in PT] , 
Here are the results: 



+ "2 Qe,nf- 

e the ! 

''"VV' + (20-2n)P'V' 



Here 'J'^^^^^p and T^^' are the symmetrizations of the tensors g^^g^'-g^^P and 



' 6 

2-3n 



2 

respectively, and J2 is the tensor 



(88-86n+n2) . . 2 (2176-768 n+82 n^W) pjk 

15 (n-4) ^ \k 15(n-4) ^ 

^-gVPj^^pkl _j_ (3n-2) (5n-54) pi,j 



2 (320-218n+27n2 



^22) 15 (n-4) y ■ fci' -r 15 

I -164-120n+45 „ij i2 2(5n-22) , k 744-250 n+31 1 ij 
+ 60 15 9 '-Ilk 15 {n-4) ^1 

2 (296-26 n+3n^) ^jfcj/ _ njklm 
15 (n— 4) 15 klm 

Here and below, for typesetting convenience, we write P^^\k^ as an alternative 
notation for V'^V^P*-' and so forth. Finally, Qg^ denotes the scalar 

(23) + iHzl^p^^.p.j _ (lizMiH) j3 + _ 6) . 

I 3n-2 I I, .« _ 8(n-6) p , i j _ 1 i j' _ 32 p .p nikjl 

We find that P%f is given by 

r 1 I lijklmpqr 

s formula, Ug 
tensors g'^i g^'' g"^P gi^ and 



In this formula, [jvklmpqr ^ji^klmp (jgj^Q|-g ^j^g symmetrizations of the 



-2n]g'^g^^g'^P - 4(n - 12)ff^J£,^'P™P, 
respectively, and 1)4 ^^' denotes the symmetrization of the tensor 

-8(-12+n) pi,- fci 4 (-12+71) (-64+5 n) piipkl , 24 n ^iy pfc; m 

5 ^ I 15 r r _4+„t/ r |m 

16(l536-530n+59n.^) 2 (480-568 n+67n^) ,,.fc; 

15(-4+n) y r- mr 15 (-4+n) » » ^rnpf 



j_±n{-6A+5n)_ ij pkl i , -96-20 n+15n^ ^ij^kl i2 , 24-5n ^t7^A:; i m 

~r 5 yrj-i- -j^Q y y -I I y y ~'\m 

8 (120-31 n+4n2) , fc; ^ 16 (l92-7n+n2) ^ ^ ^fcmZp _ 16^i j^kmlp 
5 (-4+n) y -^1 15 {-4+n) f ''^mp'-' i^'^ m p^ 



'154/ ^ mpq^ 
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-4 (-19200+8468 n-980n2+27n3) ^^,.^ ,^1 ^ 

45 (-6+n) {-4+n) f ^kl\mr 

4 (2592-5262 n+625 n2+14 1 



,3^ 



315 {-4+n) y ^f''-^ K 

4 (74928-21908 «-968 ^^+283 n^) ^ ; 

315 (-6+n) {-4+n) ^ ^^'1^^ 

8 f311616-146460n+27484n2-2167n3+5n4+7n5) ■• , ^, 

^ 315 (-6+n) {-4+n) » ^ kl^ rn^ 

2 (-44928-10224n+14400n2-3206n3+203n'*) 

45 (-6+n) (-4+n) ^kir J 

I -2560+1568 n+420n^- 105 n3 7 |3 , 14820-3650 n+231 n^ „i7 | | k 
-I- 210 y -J i- 315 y -JlfcJ 

4(-6-31n+5n2) ■ 

2 (480384-238464 n+44040n2 -3346 n3+91n4) .• 3(_48+7„) j- 
+ 3i5(_6+n){-4+n) 9'■'^^fc^J| 35 9 ' ^\k I 

I 8 (188+21 n) ^j^' p rikmlp 
315 " ~ kl\mp'^ 

8 (79488-20160 n+1020 n^+50 n3+7 n^) , r^fcmip , 16«ii r^^^/p 

315 (-6+n) (-4+n) f ^kl^mp^ 45^ ->^klmp^ 

-^g'^Cklmp\qC^'^^''{ + 16(-4+n) ^,,p^^^fc^^^^,p^g 

— iGo^ir*,, (7^^ r'^iJPf' _ 16oii/7, , nkm nlqpr 
45y ^ klmp^ q r'~' 45" ^ klmp^ qr'^ 

Figure 2. The tensor E*-^ 

We let denote the symmetrization of the tensor E*-' + F*-^ + G*-^, where 
E'^, , and G*-' are as given in Figures [21 131 and@I Finally, Ql^ denotes 
the scalar given in Figure [HI 



2.3. Branson's Q-curvature. We have used to indicate a conformally 
invariant operator between densities, '■ S[k — n/2] — > £[—k — n/2]. Sup- 
pose we choose a metric from the conformal class. Then we can trivialise 
these density bundles, and so P2fc gives an operator between functions 
on the Riemannian (or pseudo-Riemannian) structure given by the choice of 
g. If we write (^^)"' for the operator given by multiplication by (^^)"' S £[w], 
then P4/ = {^9)k+n/2p^^^^g^k-n/2kf_ rpj^g Qj^g Operators as discussed, 

for example, in are the PlJ^. In this form the operators are not invariant 
but rather covariant (see below), and conformally invariant operators are 
often discussed entirely in this setting. For many purposes the difference 
between and is rather small. In particular, since the Levi-Civita 
connection corresponding to g annihilates the formulae above for the 
also serve as formulae for the operators P^fc- From Theorem 12.51 and the 
formulae , © , lfTT|l and (|T^ , there is a universal expression for which 
is polynomial in g, g~^, C, P, V, and the V covariant derivatives of C and 
P, and the coefficients in this universal expression are real rational functions 
of the dimension n which are regular for all odd n and all n>2k. 

Let Q2k,n be the local invariant -Pf^l on a conformal n-manifold. Since 
is formally self-adjoint (FSA) it is clear we can write it in the form 
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8 (63648-43740 n+8300ri2_405n3+37n'*) ^j 
315(-6+n)(-4+n) ^^ir | 

2 (-1728+1164 n-292n2+n3) jj- ; 8 (-11108+5601 ri-447 ri2+i4 jpikl 
105 {-6+n) (-4+ji) P l'^ ' ^ 315 (-4+n) ^^l-l ^ I 

16 (-57312+18900n-1240ri2_5„3^27n4 ) pik jl 

315 (-6+n) (-4+n) ^^'^ I 

4 (-23040+98136 n-18188n2+134n3 +93 ) j ^ 

315 (-6+n) (-4+n) 
4 (-7754-1224 n+29n^) ^ pjfc « 

315 ^ ^=1'^ 

8 (589824-282624 n+54804n2-4472n3+17n'l+6n^) pifcpjZ 
^ 105 (-6+n) (-4+n) ^kir r 

4 (-51384+36938 n-9529n2+662n3+93n4 ) • ^ 

315 (-6+n) (-4+n) '=1''^ I 

4 (47232-112776 n+38488n2-4444n3+165n4) pijpkl 
"I 105 (-6+n) (-4+n) ^'^^'^ 

2 (-53720+22678 n-3194n2+63n3) ■ ^. 

315 (-4+n) ^kl\ r | 

4(l84896-132840n+30900n2-2390n3+59n4) ^ ^i., 
^ 315 (-6+n) (-4+n) "^^^'^ I 

2 (-2880-37008 n+17564n2 -2432 n3+21n* ) ^ 

315 (-6+n) (-4+n) ^ 1^= 

4 (-33408-425736 n+228064 n2-37810 n3+2094 n*+21 n^) • -g, 
315 (-6+n) (-4+n) ^ 

I -9216-1264 n+1568n^-105n3 pi, ,2 , 9336-8206 n+875 

105 r J -I- 3-^5 J| J| 

8 (-45432+19134 n-2965 n2+138 n^) j ^- ^ 

315 (-6+n) ^ 'si -"l 

2 (24768-1404 n-2584n2+215n3) -^. ^ 4(6108-1487n+70n2) ■■ ^ 

315 (-6+n) I's-'l "I 315 -^1^ 

2 (-40320-77904 n+52144n2-9826n3+651n'*) 

315 (-6+n) (-4+n) "'"'l 
4 (934848-481176 n+106800n2-11644n3+457n'') , ,1. 
^ 315 (-6+n) (-4+n) ^^-'l 

4 (-15984+8490 n-1411n2+85n=') j^- 
105 (-6+n) (-4+n) "'l ^ 



Figure 3. The tensor F'-' 



= P|^^ + Q2fc,n; where P^f} has the form ^S^j^d with 5 the formal adjoint 
of d and an order 2k — 2 differential operator. 

By setting u> = in the formulae ^ and (jHl we see that, as an operator 
on £, Da factors through the exterior derivative d. At least this is true given 
a choice of metric g from the conformal class. Thus in dimension uq = 2k it 
is clear from Theorem 12.51 that Pno is also a composition with d. Thus the 
term Q2k,n vanishes in dimension 2k = tiq. Using this and a careful use of 
classical invariant theory one can conclude that in fact Q2fc,n = Q^fc n ■ 
In the previous section we gave explicit formulae for Q^^^Q^^ Qfn- 

Clearly „ is also given by a formula rational in n and regular at 
n = riQ := 2k. In dimension no, Quq ■= Qno,no is by definition (modulo a 
sign (—1)'^) Branson's Q-curvature, and for compact conformal reo-manifolds. 
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8 (-3984+4492 n^296n^+13n3) ^ mrikjl 
315 (-6+n) (-4+n) ^kl\m ^ 

4 (33984-21624 n+4540n2 -346 n3+21n4 ) ,^ikjl 

105{-6+n){-4+?i) r'fci-J^-^ 
8 (-19392+14302 n-2403n2+128n3) ^ikjl 
315 (-6+n) (-4+n) -^1^'^ 

, 4 (-33408+13680 n-1082n2+95n3) ^ ^ifejm / _ 80^^ rikjm Ip 
315 (-6+n) (-4+n) ''^kl\m'^ | q^'^ klmp^ \ 

_ 16 (-52992+13248 n-472n2-20n3+3n'') ■ ^iklm 

63 (-6+n) (-4+n) ^kir" m'^ 

, 4 (69168-25940 n+1472n2+15n3) ^ iklm j 

315 (-6+n) (-4+n) ^kl\m^ \ 

, 8(l7088-16264n+652n^+94n3+5n4) ^ ^ i 128 (31+2 n) n i riklm 
"I 315 (-6+n) (-4+n) ^kl^ H 315 ^kl\ m^'' 

-mC'klm\pC^^^''r + ^6(-J+3n) j^i^^^^j7fcm _ 

32 (231+2 n) pi ^,mfc; , 8 (SOlli- 1570 » -oOli »-^+21 ^ k^r'nlp 

-r r fciim"-^ -I- 315 (-6+n) (-4+n) f^fcZ'-'mpL. 

,64r rikl r^mva 16 (22200-7798 n+215n^ +8 n^) ^ kr3Pl^ 
-T^^'^klmp^ 315 (-6+n) (-4+n) ^kl^ m p'^ 

+ ^CklmpC"'"'gC^''"' - MP'^CklmpC"^'" 

16 (2692-143 n+25n2) • j ^kmlp _ 16 (-33+20 n) p» ^kmlp 

315(_4+n) ^kl'-^ m p<-^ 315 ^ k'^ Imp'^ 

— 32/7 jknlmpq 



Figure 4. The tensor G*-' 



Jj^Qno is a global conformal invariant. To see this, observe that the confor- 
mal invariance of P2fc is equivalent to the covariance law 

n-\-2k ^ rt n — lh 

^-^Plk = Plk^-^^ 
where g = and we regard the powers of as multiplication operators. 
Applying both sides to the constant function 1 we obtain ^~2^ ^~^Q%^ = 

Oh y~^n „ n — 2fc ^ ^ n — 2k 

— 2 — Q2k^ ^ + ^ ^ • Expanding this out yields a universal transfor- 
mation formula for Q2k- Since is a composition with d it is clear that 
we can divide this formula by Then in dimension no = 2k we obtain 



where T = logri. Then, if we denote by Eg the volume density associated 
with a metric g, we have €g = ^'"°€g, and the conformal invariance of fj^ Q^^ 
is clear. Recall that a choice of metric g determines a canonical section of 
£^[1] by (C^)~^sr = g- It is convenient to redefine Q^f^ to be (^9)^"" times Qgfc 
as above. Then Qgfc i^ valued in n] and the transformation law simplifies 
to 

(24) Qf, = Ql, + <55f,dT, 

where 5 and are now also density valued. Note that we can also write 
this as = + P2k^ as P|^^ agrees with P^^. in dimension 2k. 
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48(-2+ra) p pij k I _ 12(-2+?i) p pjj fe i 192 n p p pifc 
-4+n I i „g+„ rtjr |A; i -I- _g^„rtjrfc/r | 

48 (-32+44 n-12n2+n3) ^ tp,-fc ; , 384(-16+3n) p pi pjfc ; 

(-6+n)(-4+n) ^tj^kl\r | "I" (_g_,_„) (_4+„) r^j r | ; r | 

I 384 (4-6 «W) i J. I I 192 (144-32 n-4n^+n3) , ^.^ 

(_6+n) (-4+n)^ ^ '^^ I' (-6+n) (-4+n)" ^'^J^'fcif' 

48 (-32+20 n-8n^+«3) ^ ^ 

H (_6+n) (-4+n) ^ij^ ^1'^ 

6 (-8+n) (256-32 n-18n2+3n3) ._ 

(-6+n ) (-4+n)" f'ii^'fcif' 
24(-176+80«-14«2+n3) ipfc/ , 4 (48+80 n-46 n^+S n^) ... 

(-6+n) (-4+n) r | -f- (_6+„) (_4+„) f^ij|fcr | J 

I 4(l2-20n+5n^) j^- fc. 384n p pi pjk, 

-6+n |fc -6+n^y^ 

96-80 n-30n^+5n^ p pij |2 , (-4+n) n (4+n) ,4 
-6+n rjjT J -t- g J 

+ (2 + 13n - 2n2) JJ|.J|^ - (-^+") (^00-38n+3n^) p^^^^.^^, 

, 24 (-2+n) (88-20 n+n^) ^ pi , fc , 16 (-48+70 «-17 n^+n^) ^ .. , 

(-6+n)(-4+n) ^ij^k\ ~l\ (_6+n) (-4+n) ""^y^ l^-^l 

I 16+4n-3n^ |2 1 i 1 -8+3 n 1 i 1 j , 2 (32-40 n+5 n^) j ■ 

_16(4-8n+n2) , j,- , 4n(-46+5n) p . . ij 

2 (-336+136n-20n2+n3) j^- 24 (IO88-6O8 n+144 n2-18 n^+n*) j 



(-4+n)" \ ' (-6+n) (-4+n)" 

+ (-38 + 5n) J| J|V - Hl_7±!!) p. .i^i^ifc + 2 (-1 + n) JJ,// 

4(-54+7n) p I ij _ I i j fc 192 (4-6 n+n^) ^ mr-ifej/ 
-6+n ^^Pl ^ -"hifc (_6+n) (-4+n)" ^■'rkl\m <^ 

I 96n p p i^ifcji 24 (-8+n) (-16-2 n+n^) 

384 (— 5+n) p p /-likjl m _ 192 p p kriiljm 
(_6+n)(-4+n)^*i^fc'l™^ I _4+n I ^ ^''"l ^ 

192 (112-16 n-6n2+n3) • ^jfcim _ _96n_ p pi A^jifcm 

(-6+n) (-4+n)" '^Jj''^fc«'"' rn"-- _4+„ f^ij | A: « |m<-^ 

I 192 p p k /~ijmlp 384 p p k /-fjplm 

' —Qj^n' ij~ kl^ m p'^ _g_|_„ r jj r^^o m 



Figure 5. The invariant 



The discussion above for Qn^ and its properties are a minor adaption 
of the arguments presented in Branson's (3- It is clear that given exphcit 
formulae for the P2k as in the previous section we can extract a formula for 
the Q-curvature as follows: Take the order part of the formula, divide by 
(n — 2k) /2 and then set n = 2k. For example Qg is obtained by setting 
n = 8 in the formula given in Figure In [SJ it is also shown that the global 
invariant is not trivial. In fact, it is established there that, on conformally 
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flat structures, Qno is given by a multiple of the Pfaflian plus a divergence, 
and so Jj^jQno is a multiple of the Euler characteristic x{^)- 

One of the keys to the importance of Q2k is the remarkable transformation 
formula We will describe a new definition and construction for Q2k and 
proof of this transformation formula. This leads to a direct formula for Q2k- 
Here to prove the transformation law we will use a dimensional continuation 
argument. (This plays a minor role and can in fact can be replaced by a direct 
proof 0). The construction is then adapted to proliferate other curvature 
quantities with transformation formulae of the general form (|2^ , and in these 
cases dimensional continuation is not used at all. See Proposition (Since 
the original writing of this Fefferman and Graham |27| have given another 
alternative construction and generalisation of the Q-curvature which involves 
the Poincare metric.) 

We work on a conformal manifold of dimension no = 2/c. For a choice 
of metric g from the conformal structure let be the section of i^^i— 1] 
defined by := (n — 2)Ya — SXa, where, recall, Ya G 1] gives the 
splitting of the tractor bundle corresponding to the metric g (as in section 
121. We can write this as a triple [I^j^g = {{n — 2), 0, —J). According 
to this definition, if ^ = fl^g then we have := (n — 2)Ya — J^a or 
= ((n — 2), 0, —J). In terms of the splitting determined by g, l\ is 
given by [l{]g = ((n - 2), -(n - 2)T„ -J - (n/2 - By © and 

= VqT this becomes = [{n — 2), — (n — 2)VaT, —J + AT), and so 

(25) li = ll-DAT. 

This observation is due to Eastwood who also pointed out (221 that on con- 
formally flat structures this yields Branson's curvature as follows. For each 
metric define Qf^ by Qfj := □2^-24. Then, by (ESJ, Q| = Qfj-n2fc-2^BT. 
Now since the structure is conformally fiat, U2k-2DA'^ = —XAP2k'^ (see 
Theorem 12.51 or e.g. [2B]). Thus we have 

g| = Qfj + XBP2fcT. 

It follows that X^Q% = X^Q% is a conformal invariant of weight nQ — 2. On 
a conformally fiat structure there are no conformal invariants of the structure 
and so this vanishes. Since this vanishes, Z^'^Q^^ is also conformally invariant 
and so must vanish. This shows that for any conformally fiat metric, = 
XbQ^ for some Riemannian invariant and also that = + P2k^ ■ 
That is, it transforms according to (|2^ . On conformally fiat structures 
one can always locally choose a metric that is fiat whence all Riemannian 
invariants vanish. Using this we deduce that is Branson's curvature, that 
is = Qlk- 

Via the theorem we can generalise Eastwood's cunning construction to the 
curved case. Note (|25|l holds on any conformal manifold. Let us define the 
operator Fc^ -.SBlk-l- n/2] £c[i- -k- n/2] by Fq^ := {k - 2)-^{n - 
2k + 2)~^D^Eck^^F)a where E is the operator defined in Theorem 12.51 
Now on a dimension tiq = 2k manifold we simply define := Fc^I^- 
From (j25p and the theorem we have immediately 



(26) 



qI = Q'+ XcP2k'^. 
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It remains to verify that for any metric g in the conformal class, Qq is 
indeed XcQgfc- In dimension n and given any metric g, if f € £[w] let 
us write DaI = D\f + wD\f, where D\f := {n + 2w -2)Za''V af - XaM 
and wDj^f is the remaining order zero part. That is, D^f = {n + 2w — 
2)YAf — XA-if- Let w = k — n/2, and assume n is odd or n < 2k. Then 
-XAP2kf = Fa^DbI = Fa^D}^/ + wFA^D^f. Let be the section 
of £"[1] corresponding to g. Recall that is parallel for the Levi-Civita 
connection of g. Since {(,^)^ is a section of £[w], we have 

XawQI,JCT = wFA^'D^iCT- 

Now X^Fa^D%{^9)w can be expressed universal expression which is 
polynomial in g, g~^, C, P, V, and the V covariant derivatives of C and 
P. The coefficients in this universal expression are real rational functions 
of the dimension n which are regular for all odd n and all n > 2k. Fur- 
thermore, from the left-hand-side of the display, this expression vanishes in 
even dimensions n < 2k and for all odd dimensions. Thus it must vanish in 
dimension no = 2k. Similarly we can conclude Z'^'^Fa^ D%{i9)^ vanishes if 
n is odd or n < 2k. Thus in dimension no, if = 0, = -D^l = 1^, 

and as Q^f, := Qn^no, we have 

Q\ = Fa^'I's = XaQI,. 
Thus we have the following. 

Proposition 2.7. Y^Fa^ 1% is a formula for Branson's Q-curvature Q2k- 

Note that if we take this as a definition for Q then the transformation 
property (|2^ arises from (|26|l which in turn is an immediate consequence of 
l(25)l . The formula itself is direct and requires no dimensional continuation. 
The only subtlety in the construction was in establishing that has the 
form X^Qgfc- We employed a dimensional continuation argument to establish 
this above, but it turns out that there is an elementary direct proof using the 
ambient construction [5]. Finally note that if we write = h^^I^, then 
//F^^4 = (n-2)Q^,. 

It is straightforward to convert this tractor formula for Q^j^ into a for- 
mula in terms of V, C, P, the metric, and its inverse. We simply expand 
Y'^Fa^I'^ using the formula for Fa^ as a partial contraction polynomial 
in Da, Wabcd, Xa, hAB, and its inverse h^^ , as obtained from Propo- 
sition 14. 5( and apply ifTTIl , and (jT^ along the same lines as the 
calculations in Section [221 In fact, as a means of checking against formulae 
or calculational errors it is prudent to calculate the entire tractor valued 
expression Fa^I^ and verify from this that only the bottom slot is not zero. 
That is that X^Fa^I% = = Z^^Fa^I^- Doing this for Q| we obtain the 
known formula 

Qi = -2PijP^^' + 2j2 - J|/. 

In terms of the Ricci curvature Rc and the scalar curvature Sc, this becomes 
= — ^Rcij-Rc*-^ -|- |Sc^ — gScjj*. For Qg the formulae are more severely 
tested by calculating Eab^^D^I^j^. For this case Eab'^^D(.I^ = + 
i^^Wa''^ DqI^ and the calculation verifies all components vanish except 
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for the coefficient of XaXb, the negative of which is 

Ql = -(8Pii|fcP'Y + 16Pi,P'V' - 32Pi,P'fcP^'= - 16Pi,P^^J+ 

8J3 - SJifc'^J + Ji/fc'^ + lePi^P^C^'^^')- 
Note that these examples agree with setting n = 4 in (^Tll and re = 6 in (EHI)- 
Using it is easy to construct examples of other functional of the metric 
that have transformation laws of the same form as ((211). We state this as a 
proposition. 

Proposition 2.8. In dimension hq = 2k, for each natural conformally in- 
variant operator Ga^ ■ 1] — > ^a[1 — ^t-o] there is a Riemannian invariant 
D^Ga^ I^B with a conformal transformation of the form 

D^GA^'li = D'^Ga'^I'b + ^T^kd'^, 

where T^^ is a Riemannian invariant differential operator such that the com- 
position 6T2j^d is a conformally invariant operator between functions and 
densities of weight —re. 

If Ga^ is formally self-adjoint, then ^T^j^d is formally self-adjoint. 

Proof. It is clear from ^ that D^Ga^I^ = D^Ga^Ib - D^Ga^DbT. 
Note that D^Ga^Db is a composition of conformally invariant operators. 
Since T is a function (i.e. is a density of weight 0), Db'^ factors through dT. 
From Proposition 12.21 it follows that the formal adjoint of D^Ga^Db also 
factors through the exterior derivative d. Thus the conformally invariant 
operator D^Ga^ Db has the form —dT^^d. 

Since ST^i^dT = —D^Ga^Db'^, the last part of the proposition is imme- 
diate from Proposition 12.21 □ 

An example in dimension 6 is to take Ga^ to be the order zero opera- 
tor \C\^6a^, where = G'''"''^Gabcd- Then D^^Ga^I^ = -4A|C|2, and 
D^Ga^I^ = -4A|Cp + 16V"|C[2V„T. We can easily make many other 
examples via the tractor objects already seen above. Other examples in di- 
mension 6 are to take Ga^ to be WacdeW^'^^^ or D^Wa^e^Df- In 
dimension 8 we could take Ga^ to be S^D^WpcdeW^^^'^ Dq and so on. 
Note all these examples have G formally self-adjoint. 

It is a trivial exercise to verify that D^Ga^Ib is always a divergence, 
and so none of the invariants from the proposition yield non-trivial global 
invariants. Thus we could adjust the definition of Q^j^ by adding such func- 
tions without affecting it as a representative of reQ^ de Rham cohomology 
and also without affecting the form of the transformation law (jSl}. Such 
changes would of course alter what we meant by S*!^, but in any case (JSl^d 
would remain an invariant operator on functions. Such potential modifica- 
tions are important from several points of view. The transformation law 
is satisfied in dimension 2 by the scalar curvature, or more precisely 
by — Sc/2. In this context it is usually called the Gauss curvature prescrip- 
tion equation. As mentioned earlier, Q2k lends itself to higher dimensional 
analogues of this curvature prescription problem. For the same reason, in 
any case where D^Ga^Ib 

is non-trivial, Q2A: + ^^Ga^Ib yields a distinct, 
and apparently equally natural, curvature prescription problem. Of course 
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then (52A; + D^Ga^Ib does not arise by Branson's construction from the 
GJMS operator P2k- But it is easily verified that it does arise via Branson's 
argument apphed to the conformally invariant operator 

P2k ■■= P2k - D^Ga'^Db : e[k - n/2] ^ £[-k - n/2], 

and according to either either construction the conformal transformation 
formula in dimension no = 2k is 

Ql, + D^GA^'li = Qlk + D^Ga^'iI + P^,.T. 

It is possible, for example, that there are settings where such natural modi- 
fications to the GJMS operators will yield operators which are positive but 
the relevant GJMS operator fails to be positive. 

3. The ambient metric construction 

The ambient metric construction of Fefferman- Graham associates to a 
conformal manifold M of signature (p, q) a pseudo-Riemannian so-called am- 
bient manifold M of signature (p + 1, g + 1). The ambient manifold M is 
Q X I, where / = (—1, 1). Henceforth we identify Q with its natural inclu- 
sion L : Q —I- M given by Q B q i-^ (^j 0) & M. Observe that Q carries a 
tautological symmetric 2-tensor go given by go = 7r*g at the point (p, g) G Q. 
This satisfies d^go = s'^go, where 6s is the natural M+-action on Q given by 
^sip^g) = {PiS^g)- We will also write 5s for natural extension of this action 
to M and denote by X the infinitesimal generator of this, i.e., for a smooth 
function / on M, Xf{q) = ■j^f{Ssq)\s=i- The metric on the ambient mani- 
fold M will be denoted h and is required to be a homogeneous extension of 
go in the sense that 

(27) i*h = go 5*sh = s'^h for s > 0. 

The idea of the Fefferman-Graham construction is to attempt to find a formal 
power series solution along Q for the Cauchy problem of an ambient metric 
h satisfying (|27jl and the condition that it be Ricci-flat, i.e. Ric(/j.) = 0. It 
turns out that only a weaker curvature condition can be satisfied in the even 
dimensional case. The main results we need are contained in Theorem 2.1 
of [26]: If n is odd then, up to a M+-equivariant diffeomorphism fixing Q, 
there is a unique power series solution for h satisfying l(27|l and Ric(/j.) = 0. 
If n is even then, up to a M+-equivariant diffeomorphism fixing Q and the 
addition of terms vanishing to order n/2, there is a unique power series 
solution for h satisfying l|27|l and such that, along Q, Ric{h) vanishes to order 
n/2 — 2 and that the tangential components of Ric(/i) vanish to order n/2 — 1. 
We should point out that we only use the existence part of the Fefferman- 
Graham construction. The uniqueness of the GJMS operators, the covariant 
derivatives of the ambient curvature and so forth are a consequence of the 
existence of tractor formulae for these objects. 

By choosing a metric g from the conformal class on M we determine a fibre 
variable on Q by writing a general point of Q in the form {p-,t^g{p)), where 
p G M and t > 0. Local coordinates on M then correspond to coordinates 
on Q. These extend (SEIISSI to coordinates {t,x'\p) on M, where p is 
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a defining function for Q and such that the curves p i— > (t, x*, p) are geodesies 
for h. In these coordinates the ambient metric takes the form 

(28) h = t'^Qijix, p)dx'dx^ + 2pdtdt + 2tdtdp. 

This form is forced to all orders in odd dimensions. In even dimensions it is 
forced up to the addition of terms vanishing to order n/2. In order, in even 
dimensions, to recover the order n GJMS operators via the procedure of [221 
we need also to assume that the metric has this form up to the addition of 
terms vanishing to order n/2 + 1. Although we only need this form to that 
order, to simplify our discussion we will assume that the form l(28|l holds 
to all orders in even dimensions too. This simply involves some choice of 
extension for the Taylor series of the components gij, and then with this 
assumption the identities discussed in the remainder of this subsection hold 
to all orders in all dimensions. We write V for the ambient Levi-Civita 
connection determined by h. 

In terms of the coordinates one has X = t-^, and if we let Q := h{X ,X), 
then Q = 2pt^ and is a defining function for Q. In terms of this we have that, 
when n is even, the ambient construction determines h up to ©(Q""/^). Let 
us use upper case abstract indices A^B, - ■ ■ for tensors on M . For example, if 
is a vector field on M, then the ambient Riemann tensor will be denoted 
Rab^ D and defined by [V^i, V^]?;'^ = Rab'^ D^^ ■ Indices will be raised 
and lowered using the ambient metric Hab and its inverse h"^^ in the usual 
way. We will soon see that this index convention is consistent with our use 
of these indices for tractor bundles. 

The homogeneity property of h in l(?7|l means that X is a conformal 
Killing vector, and in particular Cxh = 2h, where C is the Lie derivative. 
It follows that V(^X5) = hAB- On the other hand, from the explicit 
coordinate form of the metric, we have that ^bQ = 2Xb, and so ^aXb 
is symmetric. Thus 

VaXb = hAB 

which, in turn, implies 

(29) X'^Rabcd = 0. 

In terms of our notation the theorem of j2Sl (mentioned above) means 
that in even dimensions the ambient Ricci curvature Rbf can be written in 
the form 

Rbf = Q"/2^'X(bKp) + 
for appropriately homogeneous ambient tensors Kp and L'^p. In fact the 
choice to extend the metric h so that it has the form restricts Ka sig- 
nificantly. Prom we have that X'^Rac vanishes to all orders. With the 
contracted Bianchi identity 2'V^Rac = VcS' (where S denotes the ambient 
Ricci scalar curvature) this implies that Ka = XaK for an ambient homo- 
geneous function K. Although it is not strictly necessary, it will simplify 
our subsequent calculations to restrict the ambient metric a little more. An 
elementary calculation verifies that we can adjust the components gij in 
so that K = 0{Q). Thus finally we have that that in even dimensions the 
metric has the form and 

(30) Rbf = Q'''''~^Lbf 
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for an appropriately homogeneous ambient tensor Lbf- (The authors are 
appreciative of discussions with A. Cap and C.R. Graham in relation to this 
point.) 

3.1. Recovering tractor calculus. Recall that a section of £[w] corre- 
sponds to a real-valued function / on Q with the homogeneity property 
f(j),s^g) = s^f{p,g), where p & M and 5 is a metric from the confor- 
mal class [g]. Let SQiw) denote the space of smooth functions on Q which 
are homogeneous of degree w in this way. We write £{11)) for the smooth 
functions on M which are similarly homogeneous, i.e. / € S^w) means 
X'^'VaI = The construction of the GJMS operators in (321 exploits 
this relationship between £[w] and £{w). We will use here the analogous 
idea at the level of tensors on M. This is developed more fully in jT2j, and 
here we just summarise the basic ideas needed presently. 

Writing 6^ for the derivative of the action 6s, let us define an equivalence 
relation on the ambient tangent bundle by Uq-^ ~ Vq^ if and only if there is s S 
R+ such that Vq^ = s~^6gUqj^. Corresponding to this we have the equivalence 
relation on M by gi ~ g2 if and only if q2 = SsQi- It is straightforward 
to verify that the space TM/ ~ is a rank n + 2 vector bundle over M/ ~. 
Sections of this bundle correspond to smooth sections 1/ : M — s- TM with the 
homogeneity property V{5sp) = s~^6'gV{p), or they could be alternatively 
characterised by their commutator with the Euler field X, [X,V] = —V. 
We will let £^{0) {£^{0)) denote the space of sections of TM {TM\q) which 
are homogeneous in this way, and we will write £"^^{w) {£q^{w)) to mean 
£^ (8" £^ (8) £{w) {£q (8) £q (8 £q{w) respectively) and so forth. (The reason 
for the weight convention will soon be obvious.) We will write £^{w) to 
mean an arbitrary tensor power of ^"^(0) (or symmetrization thereof and 
so forth) tensored with £{w) and we will say sections of £^{w) are tensors 
homogeneous of weight w. (We use the term "weight" here to distinguish 
from the homogeneity "degree" fT2| as exposed by the Lie derivative along 
the field X .) Of course this construction is formal at the same order as the 
construction of M, but upon restriction to Q, TM / ~ yields a genuine rank 
n + 2 vector bundle over M = Q/ ~ that will be denoted by T or T"^. 

It is immediate from the homogeneity property of h that if U and V are 
sections of o (0), then the function HabU^V^ is in ^(0). Restricting to Q 
we see that HabU^V^ descends to a function on M. Prom the bilinearity 
and signature of h it follows that h descends to give a signature (p+ 1, g + 1) 
metric on the bundle T . We can use this to raise and lower indices in 
the usual way. 

Observe that G ^^(1). Thus if G ^(-1), then (/jX^ G ^^(0). 
The same is true upon restriction to Q, so we have a canonical inclusion 
1] "-^T with image denoted by T^. We write for the natural section 
of := (8 ^[1] giving this map, and so on Q, X^ is the homoge- 

neous section representing X^. Clearly then 1— > h'^^X^V^ determines 
a canonical homomorphism T and we let denote the kernel. 

Recall that Q was defined to be HabX^X^ and that this was a defining 
function for Q. Thus X^ is a null vector for the metric , and it follows 
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immediately that C T^. There is a simple geometric interpretation of 
and . Observe that T^[l] corresponds to sections of £q{1) that are 
annihilated by contraction with Xa- On Q we have that Xa = ^VaQ, 
so along Q the sections of £q{1) corresponding to T^[l] are precisely those 
taking values in TQ C TM\q and which are invariant under the action of 6'^. 
Then, since X is the Euler vector field, it follows that T^[l] corresponds to 
functions in £q{1) taking values in the vertical subbundle of TQ. Of course 
the map Q ^ M is a submersion, and so T''[l]/T^[l] is naturally isomorphic 
to = TM. Tensoring by £[-1] we have T°/T^ ^ '£'"[-1], and we can 
summarise the filtration of T by the composition series 

T = £[1]^£''[-1]^£[-1]. 

It is now straightforward to observe that the ambient Levi-Civita connec- 
tion V also descends to give a connection on T. First, from the defining 
property that V preserves the metric it follows that if G £^{w) and 
G £"^{10'), then U^VaV^ G £^{w + w' - 1). Then since V is torsion 
free, we have that ^ xU — 'VjjX — [X, f7] = for any tangent vector field U. 
So if [/ G £^{0), then ^xU = 0, as, in that case, [X, U] = —U. So sections 
of ^^(0) may be characterised as those which are covariantly parallel along 
the vertical Euler vector field. These two results imply that V determines a 
connection on T. For U ^ T let U he the corresponding section of <5g(0). 
Similarly a tangent vector field F on M has a lift to a field V G ^"^(1), on Q, 
which is everywhere tangent to Q. This is unique up to adding fX, where 
/ G ^(0). We extend U and V homogeneously to fields on M. Then we can 
form ^yU. This is clearly independent of the extensions. Since ^ xU = 0, 
it is also independent of the choice of F as a lift of V . Finally, it is a section 
of ^"^(0) and so determines a section VyC/ of T which only depends on U 
and V . It is easily verified that this defines a covariant derivative on T . 

Let us summarise. By the above construction the ambient manifold and 
metric construction of Fefferman and Graham naturally determines a rank 
(n+2) vector bundle T on M . This vector bundle comes equipped with a sig- 
nature metric /i'^, a connection V'^, and a filtration determined 
by a canonical section Xq- of T[l\. Furthermore if is a smooth tangent 
field on M and (/9 is a smooth section of <?[!], one easily verifies from the 
above that the image of v°-V'l^{ipX^) lies in and that composing with the 
map to the quotient /T^ recovers This is a non-degeneracy property 
of the connection. This with the fact that V"^ preserves the metric means 
that T is a tractor bundle with a tractor connection in the sense of 
Since V is Ricci flat it follows that V'^ satisfies the curvature normalisa- 
tion condition described in ^HIE]- (This is shown explicitly in [Ej.) From 
this and the non-degeneracy we can conclude that T"^ and are a normal 
tractor bundle and connection corresponding to the defining representation 
of SO(p + l,g + 1). That is we can take, = £^, X^ = X^, and 
to be the usual tractor connection as in Section |21 We henceforth drop the 
notation T. 

We can also recover the operators introduced in the tractor setting. Ob- 
serve that the operator Dap '■= 2XypW a] annihilates the function Q on 
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M. Thus Dap gives an operator £q{w) ^^p] ^ ^q('"^)) it is a triv- 
ial matter to show that this descends to Dap '■ £^[w] —)■ £[ap] ® £^[w] as 
defined in Section [21 (Here, of course, is the weight w tractor bun- 

dle corresponding to £q{w).) Now we can formally follow the construction 
of Da- First one calculates that, for V G £^{w), and using (|29|) . we have 
h^''DAiQD\B\P),V = -X^gDpyy, where 

(31) DaV = {n + 2w -2)VaV - Xa^V, A := V^Vs- 

Then we observe the map £p{w — 1) ^ £(^pQ-)g{w) given by Sp i— X(Q5p)g 
is injective. It follows immediately that, along Q, pT]! is determined by 
the equation h^^ Da(qD^^^p-^^V = —X (^gOp^^V and so is precisely the 
operator Da '■ £q{'w) — > £^ (^g(t« — 1), which descends to Da '■ £^[w] — > 
£a ® £^[w — I]. In particular this is true when w = 1 — n/2, and so A : 
£^{1 — n/2) — £^{—l — n/2) descends to the generalised Yamabe operator 
□ : £■*[!- n/2] ^ £:*[-!- n/2]. We will take (EU as the definition of 
Da on M. Although we will not need it here, let us point out that Dap as 
defined above acts more generally on sections of tensor bundles on M and 
not just sections which are homogeneous. Following through the argument 
above in this more general setting yields a generalisation of the operator Da 
on tensor bundles given by Da = nV a + 2X^'W p^ A — Xa^- This still 
has the property that along Q it acts tangentially. 

Observe that h^^ D a[qD\p\p-^^V is only of the form —X(^QDp-^^V to 
order along Q and that although Da acts tangentially to Q to this order, 
it does not commute with Q. In fact for any tensor field F, homogeneous of 
weight If on M, from l(3T1) we have 

(32) DaQV = QDaV + AQVaV. 

So, along the Q = surface Q, Da acts tangentially, but. Da does not act 
tangentially to other Q = constant surfaces. Nevertheless this allows us to 
conclude that if U and V are tensors of the same rank (and with U + QV 
homogeneous of some weight), then 

Da, Da,{U + QV) = {Da, ■■■Da,U) + QW 

for some tensor W. Thus, along Q, Da, • • • DaJJ is independent of how U 
is extended off Q. The identities 

(33) XaD^V = w{n + 2w-2)V -Q^V 
and 

(34) D^XaV = {n + 2w + 2){n + w)V -Q/^V 

will also be useful. Here is a tensor which is homogeneous of weight w. 

We are now in a position to show directly how the tractor field Wabcd 
is represented in the ambient setting. Let us for the while restrict to n / 4. 
Note that the curvature of the ambient connection Rabcd is a section of 
^^a_bcd(— 2) and so determines a section of the tractor bundle £abcd[—'^]- 
We will write Rabcd to denote this section. Let V G £^{w). From we 
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obtain 

[Da,Db\V = {n + 2w -2){n + 2w -4.)[Va.^b\V 
-2{n ■^2w- 2)X[^[A, V s^Y . 

Now let G E^. We write V = G "^2(0) for the corresponding field on 
Q, and extend this homogeneously to a field on M. Then, along Q, we have 
(see remark below) 

\Da, Db]V^ = {n-2){n- A)Rab^eV^ + 4(n - 2)X [aRb]f^ e'^^V^ ■ 

Thus, since Rbfce = = X p'V^V^ , this implies 

[Da, Db]V^ = (n - 2)(n - 4)Rab^eV^ + 4(n - 2)X[aRb]f^ eZ^ jV^V^ . 

Comparing this with (|T2|l (with w set to in that expression) we can at 
once conclude that X[aWbc]deV^ = (n — 4)X[aRbc]deV^ ■ Since this 
holds for any section of it follows from the definition of Wabcd that 
X\^A^BC]DE = ^[aRbc]de- Contracting with f we have immediately 
^[aRb]f'~^ eZ^ f = X^A^B]F^ eZ^ f. Substituting this in the above display 
and once again comparing to (|T2|l we now have that WbcdeV^ = (n — 
4:)RbcdeV^ for all V^, and so 

(35) Wbcde = (n - 4:)Rbcde- 
Remark: Note that 

[A, Vb]Vc = 2RebcfV^V^ + {V^Rebcf)V^ + Rbf'^^Vc. 

From the contracted Bianchi identity 'V^ Rebcf = 2'V^cRf]b, so in odd 
dimensions the last two terms of the display vanish to all orders. In even 
dimensions recall we have that, along Q, Rbf vanishes to order n/2 — 1 and 
so in all even dimensions, other than 4, these last terms also vanish along Q. 

4. The GJMS operators 

Using the properties of Da, we observed in the previous section that if V is 
a tensor homogeneous of weight 1 — n/2 then, along Q, AF is independent of 
how V is extended off Q. So A gives an operator A : ^*(l-n/2) ^ ^*(-l- 
n/2), and this descends to the generalised conformally invariant Laplacian 
(or Yamabe operator) as in (CU)- The observation that the conformally 
invariant Laplacian (on densities) can be obtained from an ambient Laplacian 
in this way goes back to ™ the conformally flat dimension 4 setting and 
to jSn] for the general curved case. For the generalised conformally invariant 
Laplacian we can also show this directly using the result 

(36) VaQ = 2Xa 

from above. Prom this it follows that if [/ is a tensor field homogeneous of 
weight w (i.e. U € £^{w)), then 

(37) [A,Q]U = 2{n + 2w + 2)U. 

Thus if F €<?*(! — n/2) and [/ is a tensor of the same rank and type but 
homogeneous of weight —1 — n/2, then 

(38) A{V + QU) = AV + QAU. 
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So clearly A.V is independent of how V extends off Q. In [53, Graham, 
Jenne, Mason, and Sparling establish a remarkable generalisation of the re- 
sult for densities which we state here in our current notation. 

Proposition 4.1. For n even and /c G {1, 2, • • • , n/2} or n odd and k G 

let f € <?q(/c — n/2), and let f € £{k — n/2) be a homogeneous extension 
of f. The restriction of A'^f to Q depends only on f and the conformal 
structure on M but not on the choice of the extension f or on any choices 
in the ambient metric. Thus there is a conformally invariant operator 

A'' : Sgik - n/2) ^ ^q(-A; - n/2), 

and this descends to a natural conformally invariant differential operator 

P2k ■■ £[k - n/2] £[-k - n/2] 

on M. 

As mentioned in the introduction, we call the operators P2k the GJMS op- 
erators. 

In this section we will describe a way that one can directly rewrite these 
operators in terms of Da, Xa, the curvature R, and just one A. As ob- 
served above, each of these corresponds to an object in the tractor calculus. 
Before we begin we need one more result from jS2]- (This follows from 
Proposition 2.2 and Section 3 from there). 

Proposition 4.2. For n even and A; G {1, 2, • • • , n/2} or n odd and k S Z_|_, 
let f G £Q{k — n/2). Then f has an extension f € £{k — n/2) uniquely 
determined modulo 0{Q^) by the requirement that Af = modulo 0{Q^~^). 
The extension is independent of any choices in the ambient metric. 

We are ready to consider an example. Let / G f (2 — n/2), and let / 
denote the section of f [2 — n/2] that it determines. Consider ADAf = 
A(2Va/ — XAAf). Since in all dimensions the ambient Ricci curvature 
vanishes along Q, we have [A,Va]/ = 0. So with the operator equality 
[A, Xa\ = 2'W A we immediately see that 

ADaJ = -XaA^I 

Thus D-D^/ = —XAPif, where P4 is the fourth-order GJMS operator (which 
agrees with the Paneitz operator). Note that according to the earlier propo- 
sition above, the right-hand side is independent of how / extends off Q. So 
the left-hand side is likewise independent of the choice of extension. In fact 
this is already clear from (p?^ and ((SHI- 

This suggests attempting to recover the higher order GJMS operators from 
AD A ■ ■ ■ Dsf ■ On conformally flat structures this is immediately successful. 

Proposition 4.3. On conformally flat structures, if f € £{k — n/2), k G Z+, 
then 

ADa,_, ■■■DaJ= {-1)'''Xa, ■ ■ ■ XA,_,A^f. 

Proof. We are only interested in local results and differential operators. So 
without loss of generality we suppose that we are in the setting of the flat 
model for which the ambient space is simply M"+^ equipped with the flat 
metric h given by a fixed bilinear form of signature {p + 1, q + 1) and the 
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standard parallel transport. The latter also gives the ambient connection in 
this setting. In the standard coordinates, X = X^d/dX^ at the point X^ , 
and the identities of the previous section hold as genuine equalities rather 
than just formally. 

We have the operator identity [A,Xa] = 2Va on sections of £^{w). 
Since the structure is conformally flat, we also have [A, ^a] = 0. It follows 
that [A"',Xa] = 2mA"'~^VA- Thus if / G + 1 - n/2), we have 

-A"' Da} = -A^[2mWA~f - XaAJ] = X^A'"+i/. 
The proposition now follows by induction on k. □ 

To relate A.DAf._-^ ■ ■ ■ Da^J and A'^/ in the general case we must take 
account of the curvature of the ambient manifold. Since this is Ricci flat 
we have that if Vb G Sa{w), then [A, Va]Vb = -2Ra^ b^'^ pVq. More 
generally if Vbc-e G £bc-e{w), then 

[A,V aWbc-e = 

(39) -2Ra^ B^^ pVqc-e - 2Ra^c'^VpVbq...e - ■ ■ ■ 

-2Ra^e'^'^pVbc-q- 

In even dimensions the ambient metric is only Ricci flat and determined 
by the conformal structure on M to finite order, as described above. For 
example for even n (Enj only holds mod 0(Q"/2-2^ 

lor mod 

has rank 0). For simplicity in the following discussion we will often ignore this 
point and assume the given calculations do not involve sufficient transverse 
derivatives of the ambient metric to encounter this problem. We will return 
to a careful count of tranverse derivatives later in the section. We will also 
henceforth restrict to n 7^ 4. This also simplifies matters. And there is no 
loss, as the results for n = 4 have been obtained above. 

It follows from the last display that if / G £{w) (and £ < n/2 if n is even), 
then 

AVAr--^Aj = 

-2i?A,%_,«VpVQ Va,_, ■■■VaJ---- 

-2VA,i?A,_/A,_,« Vp Vq Va,_3 ■■■VaJ---- 

-2Va, ■ ■ ■ VA;RA/A,'^VpVQf+ 

Va,---Va,A/, 

where here all V^'s act on all tensors to their right and the result is mod 
0(Q"/^~^) if n is even. We may apply the Leibniz rule to (|4n|l . The term 
V AiRai_^^ Ae_2^'V p'V qV At_3 - - - Vai/, for example, becomes 

(Va,-Ra,_/a,_2'3)VpVqVa,_3 • • • Vai/+ 
i?A,_/A,_2^VA,VpVQVA,_3 • • • VaJ. 
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Often we will not require the details of contractions or the value of coeffi- 
cients, and so we might write the last result symbolically as V/?V^~^/ = 
(Vi?)V^"V+-RVV- (In this informal notation we will write V to indicate 
a Va which is not part of a A. For example, it may have a free index or 
be contracted to the ambient curvature R.) We may repeatedly apply the 
Leibniz rule to l|irHl in this way until all of the terms on the right-hand side 
are of the form (omitting indices) (V^i?) V^/. We might write the result 
symbolically as 

(41) AVV~ = V^Af + ^iVPR)^/. 

Note that each term of the second sort on the right-hand side has q >2 and 
p + q = i. Although in these symbolic formulae we omit the details of the 
contractions and the coefficients, we really want to regard these expressions 
as representing precise formulae. The idea of this notation is simply to 
manifest explicitly only the aspects of the formulae that we need for our 
general discussion. 
Now observe that 

^n + 2w-2i- 2)Va,+i Va, • • • VaJ" = 

Da,+, Va, • • • / + Xa.^^AVa, • • • /, 

or, in our symbolic notation, {n + 2w-2£- 2)V^+V = DV'^f + XAV^f. 
We can substitute into the right-hand side of this and so observe that 
if n + 2w — 2£ — 2 ^ 0, then we can replace a term V^"*"^/ by the expression 
i:>VV+^E(V^i?)V''/ + ^V^A/. Supposes; = A; -n/2. Thenn + 2u;- 
2£-2 = 2{k-e- 1), and we have 

(42) 2{k-£- 1) V^+V~ = DV^f + X ^( VPfl) VV + X V^A/. 

In each term of the sum we again have q > 2 and p + q = i. Note that the 
left-hand side of has at most i + 1 transverse derivatives of /. Apart 
from the term XV^A/, which we will deal with below, the right-hand side 
has at most i transverse derivatives of /, as D acts tangentially to Q. Our 
strategy below will be to replace V's with D's beginning from the left. 

We may apply similar reasoning to R. Since R has weight —2, we have 
{n-2m-4:)V"'R = DV"'-^R+ X AV""'^ R. Here we have used the same 
informal notation that we used with /, above. By (|H9|1 we may write this 
as (n-2m-4)V'"i^ = £>V™-ii? + X^(VPi?)V5i? + XV™-iAi?. Now 
note that since R is Ricci flat, we have 

(43) A ^J^BCDE = 

2 {R cb^ Rfade + R CD^ Rbafe + R ce^ Rbadf) , 

from the Bianchi identity. In odd dimensions this holds to all orders. In 
general we have ARbcde = 2(VBV[£)i?£;]c- VcVpi?^]B)-FO(i?^) where 
0(i?^) indicates the quadratic term in the display. Using, once again, that 
in even dimensions Rab = Q"'^'^~^Lab it follows that (V b^[d^e]c ~ 
^C^ID^^eib) vanishes to order n/2 — 3, and so holds to that order. 
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Thus we get the simphfication 

(44) (n - 2m - 4) V^i? = DV'^-^R + X ^(V^l?) V«i^, 

where in each term of the sum, p + q = m — 1. In even dimensions we need 
m < n/2 — 2. This follows immediately from the previous paragraph. That 
is, we need (n — 2m — 4) > 0. 

Our effort to relate ADAf,_^ ■ ■ ■ Da^J and A'^/ involves another identity, 
viz 

(45) A(V*A"i?)S = (AV*A"i?)S + (V*A"i?)AE + 2(V*+^A"i?) VS. 

Here E is any expression (for a linear operator) which, in terms of our in- 
formal symbolic notation, is a polynomial in V, A, R, and /. We also need 
the following fact which follows from the above: 

Lemma 4.4. Suppose n is odd or t + u < n/2 — 3. Then on Q there is an 
expression for V* A"i? as a partial contraction polynomial in Da, Rabcd, 
Xa, hAB, and its inverse h^^ . This expression is rational in n, and each 
term is of degree at least 1 in Rabcd- 

Proof. Repeatedly use (EH), HSI), and (g^ to rewrite V*A"i? as a sum of 
terms of the form (V^i?) • • • {W"^ R). In doing this we convert some A's 
into pairs of V's via but at most one V from each pair acts on any 
given R. Thus in even dimensions, Vi < n/2 — 3, i € {1, • • • , j}, and we may 
construct the desired partial contraction polynomial by repeatedly applying 
(|331l to the terms (V^ii) • • • {V^^R). In even dimensions, using (jtZl and 
Va<3 = 2Xa with the restriction t + u < n/2 — 3, we see that that (jSHI); 
iHSjl, and (|331l all hold to sufficient order. □ 

Now let / G £Q{k — n/2). Suppose / € £{k — n/2) is any homogeneous 
extension of / as in Proposition 14.21 We will consider A.Da^_-^ ■■•DAif, 
where k is a positive integer. If n is even, we assume that k < n/2. Let 
us systematically rewrite this in terms of (— 1)'^~^X^^ •••X^^^A^/ and 
curvature coupled terms via the following steps: 
Step 1: Observe that 

ADa,_, ■■■DaJ= A(2V^,_, - A) . • • {2{k - 1)Va, - X^, A)/. 

Expand this out via the distributive law without changing the order of any 
of the operators. 

Step 2: Move all X's to the left of any V or A via the identities [V^, X b\ = 
hAB and [A, X^] = 2'V a (which hold to all orders). 

Step 3: Move all A's to the right of any V's (other than those implicit in 
A) via and In even dimensions one of course needs to be careful, 
since l|in|l is vahd only if £ < n/2 and holds mod 0((5"/^~^). Elementary 
counting arguments (along similar lines to the discussion in the next para- 
graph) quickly establish that for terms encountered we have i < k satisfied 
and with no more than {k — £ — 1) transverse derivatives of the result. Since 
we assume k < n/2 when n is even the use of (|ir)|l is valid. Next by the 
proof of Proposition 14. 3( we may cancel all terms not explicitly involving 
the curvature except for the term {—1)''^^Xai • • • ^A^.i A'^/. (The proof 
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of Proposition EHl involves only the identities used in Steps 1 and 2 with just 
the difference that these are applied in a different order.) We thus obtain 

(46) {-l)^-^X''-^A''f + h'X''{VP'A''^R) ■ • • (VP'* A*"-^!?) VA^^, 
where d > 1 in each term of the right-hand part. 

At this point let us take stock of what we have. For each term in the result 
of Step 1, the sum of the number of A's in the term and the number of V's in 
the term is exactly k. In Steps 2 and 3 some A's may have been exchanged 
for V's via the identity [A,Xa] = 2Va or for R's via the commutator 
[A, Vyi], and similarly we may have lost some V's by [Vyi, X^] = Ixab- On 
the other hand, we may have converted some A's into pairs of V's via 
note that at most one V from each pair acts on /, and similarly at most one 
V from each pair acts on any given R. Thus for each term of the sum in 
we must have d + q + r < k. Since d > 1, it follows that k — q — r > 1. 
Note that each R in (|in|l is followed by at least two V's. Thus at each step 
in the construction of the right-hand part of iP6|l , each of the rightmost two 
V's of each term arose from Steps 1 and 2, and not from the use of (|i5|l . 
It follows that at least two of the V's in V^A*^/ did not arise from (|i5|l . 
Thus g > 2, and for any i € {1, • • • , d}, Pi + ri + 3 < k. Now suppose n 
is even. Then, by assumption, k < n/2, and for i € {1, • • • , d} we have 
Pi + fi < n/2 — 3. Since the ambient metric is determined modulo terms of 
0(Q"'/^), it follows immediately that the metric connection V is determined 
modulo terms of 0((5"/^^^). Its curvature R is similarly determined modulo 
0((5"'/^~^). Now when A = V^Va acts on functions, its rightmost V is 
really just the exterior derivative. Thus as an operator on functions, A is 
determined modulo terms of 0(Q"/^~^). It now follows from ^ and ^ 
that, as an operator on f (A;— n/2), all terms in the sum of are determined 
uniquely modulo 0(Q). If n is odd, the ambient metric is determined to 
infinite order so certainly the same is true. 

Next, by Proposition 14.21 we can assume that / satisfies A/ = modulo 
0{Q^~^), and given /, this determines / uniquely modulo 0{Q^). This will 
simplify our arguments. The end result will be independent of this choice. 
Since A; — g — r>l, we see immediately that all terms in (|i6|l with r > 1 will 
vanish modulo 0{Q). We will thus delete these terms. From the inequality 
k — q — r > 1 and, in even dimensions, the inequality Pi + ri < n/2 — 3, it 
follows that we can carry out the next step. 

Step 4: First rewrite as 

(47) {-l)''-^X''-^A''f + /i"X^(VPi A'^iil) • • • (V'' A''rf/?)VV- 

Then repeatedly use and Lemma to eliminate all V's and A's from 
the right-hand part of this expression. The use of introduces additional 
A's. But terms containing these A's vanish modulo 0(Q), and we cancel 
them as soon as they appear. We obtain as result 

(48) {-l)^~^X^-^A^f = AD^-^f + Y'^'^f^ 

where, in terms of our informal symbolic notation, the operator ^' is a poly- 
nomial in X, D, and R. The exponent s here is not claimed to bear any 
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relationship to the s from earher. The only differential operator of non-zero 
order used in the formula is D. Thus although we used the / satisfying 
A/ = modulo 0(Q'^~^) to obtain (@Hll, observe now that it follows im- 
mediately from l(32|l that each term depends only on / and is otherwise 
independent of the extension /. Thus for any extension /, (|i8|l holds mod- 
ulo 0(g). 

Remark: At this point it is worthwhile to justify our use of in Step 4. 
First note that in each term in l|T7|l we have q < k — 1, hy the counting given 
above. Thus in £ + 1 will always be at most k — 1, i will be at most 
k — 2, and k — £ — 1 will be nonzero. We may therefore solve for V^"*"^/ in 
On the other hand the use of may generate additional curvature 
terms (V^i?) V^/. But p + q = £, where q > 2. Thus in even dimensions, 
p<£-2<k-4< n/2 - 4, and we may apply Lemma lO to V^R. 

In the final step we will use the fact that (n — 4)i? descends to the tractor 
field W, X descends to X, h descends to h, and that A : £"*(! — n/2) — > 
S^{-1 - n/2) descends to □ : f *[1 - n/2] £■*[-! - n/2]. 
Step 5: In the right-hand side of (|i8|) make the following formal replace- 
ments: / with /, A with □, X with X, h with h, R with W/{n - 4) 
(in dimensions n ^ 4) and D with D. The result is a tractor formula for 
(_l)fc-iX'^~ip2fc/- We state this as a proposition. 

Proposition 4.5. There is a tractor calculus expression for the GJMS op- 
erators of the form 

(49) XA,---XA,_,P2kf = 

(-l)^-inZ)A,_, ■■■DaJ + ^A,^,-A,'''^DpDQf, 

where f £ £[k — n/2] and ^ is a linear differential operator 

: £pQ[k - 2 - n/2] ^ £a,_,-A, [-1 - n/2], 

expressed as a partial contraction polynomial in Da, Wabcd, Xa, It-ab, and 
its inverse h^^ . The expression for ^ is rational in n, and each term is of 
degree at least 1 in Wabcd- 

Proof. It is clear from the argument of this section that 

Xa,--- XA,_,P2kf = {-l)'-'aDA,_, ---DaJ + ^a,.,-aJ, 

where "i^Ak-i- Ai is a linear differential operator on / expressed as a partial 
contraction polynomial in Da, Wabcd, Xa, hAB, and its inverse h^^. It 
is also clear that that this expression for is rational in n and that each 
term is of degree at least 1 in Wabcd- Furthermore, recall that in Step 4 
we used (|i2|l and Lemma ^3 to convert the expression V^/ of (|i7|l into an 
expression in D, X, R, h, and . Since q > 2 in (|i2|l and (|i7|l . it follows 
that each term of this tractor expression ends in two consecutive D's. The 
result now follows. □ 



We conclude this section with examples. 
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4.1. Examples. The simplest example of our procedure is the Paneitz op- 
erator -P4, which we treated at the outset of this section. Recall that we 
obtained ODaI = —XaPaJ, and it is clear that the tractor expression on 
the left-hand side of this is independent of any choices in the ambient con- 
struction. This is as guaranteed by the argument following Step 3. 

The next simplest case is of course the operator Pq. By assumption then, 
n / 4. Let / denote a section of£^[3 — n/2]. Let / be a section of £^(3 — n/2) 
such that its restriction to Q agrees with / and such that, A/ = Q^g for 
some smooth g G £"(— 3 — n/2). Expanding out ^DaDbI according to 
Steps 1 and 2 gives 

/^DaDbJ = 

XaXbA3/ + 2Xb[V^,A]A/+ 

2Xa[Vb, A] a/ + 4XaA[Vb, A] / - 8[V^, A] V^/. 

Since [A,V^] vanishes on functions mod O(Q^), the third step reduces to 
ADaDbI = XaXb^^J - 8[Va, A] Vb/ = 

along Q. The fourth step is simply the observation that on Q (with / as 
above) we have 

SVcVeI = DcDeI 

and so 

XaXbA^I = ADaDbJ + 2RA^B''DcDEf. 

As we have observed generally, at this stage none of the terms on either side 
depend on how / extends off Q. Thus finally we have 

UDaDbJ + -^WA^^B^DcDEf = XAXsPef, 
n — 4 

where Pq is the sixth-order GJMS operator. Thus as promised, we have 
recovered the tractor formula found by other means in Section O 

Our final example is Pg. By following Steps 1 through 4, above, and by 
applying we obtain 

-XaXbXc^^I = 

ADADBDcf + 2RA''B'^DpDQDcf + 2RA^c^DpDBDQf 
-j^^XA{DERB''c'^)D''DpDQ~f 

+AXARB^C^Rp^Q^DEDFf - 2XAUB^C^DpDQf 

~^XAXEUB''c^D''DpDQf 

+ j^^XAX^{DERB''c'^)Rp^Q''DFDGf. 

Here U b^c^ denotes the tractor field 

2 {R^^b^Reac'^ + R^'^'c'^Rbaf'^ + R^'-'^^Rbacf) • 
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To demonstrate explicitly that Ps is formally self-adjoint, a variation on 
this formula is preferred. It is a straightforward exercise to rewrite the above 
equation as follows. 

-ADaDbDcJ - 2RA''B'^DpDQDcf - 2RA''c'^DpDBDQf 
-^XA{UB''c^)DpDQf + ^XAD''RB''c^DEDpDQf 
+4^XARB''c'^Rp^Q'^DsDTf. 

This together with yields the tractor formula of Proposition 12.31 
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